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INTRODUCTION 

BEACON  denotes  a  set  of  programs  organized  to  form  a 
computing  tool  to  study  defense  of  the  continental  United  States  against 
strategic  nuclear  attapk.  The  work  has  been  predominately  supported 
by  the  USAF,  Headquarters  Systems  Command  under  Contract 
AF  18(600) -2793.  The  purpose  is  not  to  evaluate  current  force  structures 
(although  this  is-hotan  unreasonable  application)  ,  but  rather  to  .obtain  a  ^ 
feel  for  the  possibilities  of  newer  technological  concepts.  Wild  fresh 
ideas  may  form  the  other  end  of  the  application  spectrum,  though  it 
is  usually  advisable  to  trim  them  down  first  by  traditional  study  pnases; 

Thus,  BEACON  has  not  been  designed  toward  the  study  of  a 
particular  weapon  system  or  even  type  of  systems,  but  (in  intent  at 
least)  toward  embracing  the  whole  set  of  potentially  interesting  intern- 
acting  systems  of  the  near  future.  Since  it  is  difficult  for  a  small  set 
of  people  to  be  good  visionaries  over  sUch  a  vast  realm,  the  method  of 
procedure  has  been  to  concentrate  bn  fundamental  logical  or  physical, 
processes  of  attack  and  defense  on  a  global  scale  and' to  keep  the 
structure  mpdular,  replaceable,  and1  open-ended,  it  has  hot  been  easy. 

A  distinction. must  be  made  between  the  general  procedure  that  may  be 
termed;  BEACON  and  any  particular  r e pr e s entatioh  of  it  that  may  exist 
at  a  particular  time  such  as  BEACON  I  dr  II  or  III,  etc*  It  ia  inevitable 
that  representations  be  incomplete  or  unsuitable  for  the  *!next” 


application.  Success  is  often  a  combination  of  existence  and  adaptability* 
In  this  light  the  representation  included  here  can  only  be  regarded  as  a 
transitory  example  of  the  general  procedure. 

BEACON  is  not.a  model,  atleast  not  according  to  the  usual 

o 

thoughts  that  the  word  model  evokes.  The  germanic  concept  that 
disingtiishes  a  model  is  the  occurrence  of  a  central  ur  main  program  i 
that  controls  the  processing  flow  according  to  some  built-in  logic  which 
imitates  or  captures  the  essence  ot  some  real  life  activity.  A  key  word 
lies  in  built-in.  The  fact  that  most  models  permit  selective  control  or  . 
parameter  variation  through  inputs  at  rim  time,  cannot  destroy  the 
essential  feature,  of  a  definable  pre  -understood  logical  path.  By  contrast 
BEACON  substitutes  for  the  central  program  a  monitor  which  controls 
during  execution  other  semantic  component  programs  for  the  user  in  the 
way  the  iiser  specifies.  The  entity  takes  the  form  of. a  computation  by 
request,  in  the  order  of  request,  and  to  the  extent  requested.  It  is  not 
farfetched  to  say  that -BEACON:  passes  on  to  the  user  the  joy  (and  the 
task)  of  "writing"  his  own  main  program  at  each  execution.  The  method 
is  simple  and  flexible  enough  to  suit  most  user  type  needs  without  the 
elaborate  procedure  a  programmer  uses  of  actually  writing  and  com¬ 
piling  a  central  program .  The  -user  pays  an, initial  -price  in  'onuer standing 
and  organizing  his  problem,  but  gains  a  degree  <St  freedom  and  parties- 
patiph.  Although  the  subject  matter  obviously  prohibits  classifying 
BEACON  as  a  general  purpose  software,  within  its  scope  this  Is  its  spirit 
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SECTION  I 

Since  BEACON  does  not  have  a  prescribed  logic,  one  cannot  v,  < 
draw  a  very  informative  flow  diagram.  It  seems  best  in  this  circum- 
s  stance  to  discuss  each  of  the  main  programs  and  then  to  illustrate 
some  ways  they  may  be  exercised  on  defense  problems.  The  function 

t 

of  the  monitor  BEAM  is  only  to  aid  iri  effecting  the  latter  and  will  be 


discussed  later. 

ATLAS 

Almost  always  ATLAS  will,  be  the  first  program  a  user  will 

wish  to  perform.  The  job  which  is  accomplished  by  ATLAS  is  to  initially 

set  up  in  core  memory  what  is  called  a  facility  record.  The  information 
’ »  • 

content  of  the  facility  record  is  to  be  a  description  of  the  status  of  the 
*  geographical  and  physical  characteristics  of  the  globally  fixed  elements 

»  essential  to  the  battle,  included  in  the  record  will  be  space  for  data 

f  from  some  of  the  other  programs  .  The  elements  are  grouped  into  sets 

with  each  set  given  a  name  by  which  it  is  referred.  What  facilities  are 
necessary  is  solely  dependent  upon  the  problem  being  studied  and  the 
programs  to  bf  used  subsequently.  Typically  there- will, be  sets  to  define 
the  targets  of  the  attacker.  These  may  be/couhterforce,  countervalue,  or 
defense  sites.  US  offense  bages  ar^  often-tncluded  fori-this  reason^,  -If—- . 


there  is  to  be  ah  active  defense  there  willbe  one  or  more  interceptor 
and/or  sensor  sets  to  affect  the  various  battles#  Still  a  third  reason 
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for  inclusion  of  a  set  may  be  for  damage  assessment  comparisons., 

ATLAS  decks  for  the  213  US  and  .26  Canadian  -urban  areas  have  been 

-  ’  ‘  v  ’  -  * 

*■  j.f 

constructed  since  they  are  so  commonly  used  for  damage  assessment. 

Thus,  tht  facility  record  Is  an  opan-«nd«d  u»sr  defined  concept. 
The  ATLAS  data  deck  required  to  supply  the  detailed  information  is  by 
far  the  largest  part  of  the  typical  input  date,.  The -common  way  of 
generating  such  a  date  deck  is  to  add >  delete,,  or  otherwise  shuffle  a 
previous  deck.  With  this  in  mind  BEACON  is  so  constructed 'that  it, is 
not  necessary  for  the  user  to  preserve  or  rememberteny  special  order  • 
nor  to  supply  counts  either  of  the  sites  within  a  set  or  of  the  sets  them¬ 
selves.  It  is  incumbent  upon  the  user  to  remember  the  names  he  .has 
used.  Another  feature  which  makes  it  easier  to  generate  an  ATLAS  deck 
from  cannibalizing  existing  decks  is  the  personal  order  and  format 
specification  each  facility  set  is  permitted.  Source  data  may  be  presented 
in  a  wide  variety  of  forms  and  more  often  than  not  the  formats  are 
mixed  in  any  given  rim.  It  is  then  possible  to  input  data  for  a  set  of 
airfields  as  :  name,  longitude,  latitude-,,  while  for  urban  areas  the 
list  may  read:  population,  name,  latitude,  longitude,  area,  .4  . 
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in  part  upon  the  nature  of  the  set.  Currently  five  distinct  kinds  of 
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facilities  are  recognized  in  BEACON  II:  sites,  blocks,  places*  barriers, 
and  meshes.  Sites  are  simple  homogeneous  elements  which  differ  only 
in  location  and  value,  such  as  may  suffice  for  many  types  of  military 
bases.  The  sites  of  a  given  set  are  assumed  to  have  the  same  physical 
‘radius  and  hardness  for  targeting  and  damage, assessment  purposes. 
Blocks  may  be  used  conveniently  when  a  cluster,  of  point,  targets  are 


close  enough  together  to  be  designated  as  a  single  squadron  or  farm  and 
yet  far  enough  apart  to  require  separate  targeting*  as  for  example  , 
hardened. missile  silos.  Each  block  is  specified  as  to  location*  size 
and  unit  value.  When  required  each  block  of  a  given  set  will  be  auto¬ 
matically  deployed  into  a  squarish  array  of  the  specified  number  of  points 
around  the  center  of  the  squadron  "headquarters".  The  hardness  and 
distance  between  points  are  set' characteristics.  ip.lacesare  used  to 
designate  those  elements,  notably  urban  areas,  which  cannot  be  reason¬ 
ably  regarded  as  the,  same  size.  The  individual  information  in  this  case 
includes  not  only  location  and  value  (population),  but  specific  area  in 
square  nautical  miles .  For  the  three  kinds  of  facilities  so  far  sites, 
blocks,  and  places  it  is  optional  but  recommended  that  analphanumberic 
designation  be  also  stored  to  identify  later  printouts. 


Barriers  and  mes 


- - 


e  imagiasi-y  ■iaciiiueoji  that  i's  they  are 


not  targetabie  nor  assessable.  Their  purpose  is  to  allow  the* user  to 
introduce  a  topology  into  the  bfttle  £n  a  form  suitable  to, his  problem. 


A  barrier  is  formed  by  simply  giving  a  list  of  fence  post  locations,  i.;?. 
latitudes  and  longitudes.  In  this  case  the  order  of  the  list  is  meaningful. 
ATLAS  will  automatically  construct  the  barrier  by  connecting  great  circle 
arcs  between  successive  posts  and  will  close  the  ring  by  connecting  the 
last  post  back  to  the  first.  As  the  barrier  is  traversed  in  ascending 
order,  ine  arcs  are  numbered  the  same  as  the  posts  being  approached  and 
outward  direction  is  taken  to  be  to  the  left.  Barriers  often  consist  of  a 
dozen  or  so  arcs  and  may  represent  such  things  as  probable  detection 
lines,  extent  of  fighter  coverage,  coastlines,  etc.  A  mesh  is  an  exten¬ 
sion  of  the  barrier  concept  by  adding  to  each  post  a  directional  heading; 
Headings  are  specified  as  degrees  clockwise  from  north  and  normally 
they  are  positioned  to  divide  the  outward  sector's  formed  by  the  adjoining 
arcs.  The  exterior  of  the  ring  is  thus  subdivided  into  sections  which 
are  numbered  the  same  as  the  corresponding  arc  (see  illustration). .  ATLAS 
constructs  great  circle  arcs  so  as  to  be  tangent  to  the  headings,  and  they 
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form  the  global  extensions  of  the  headings..  ;(A  heading  of  due  EAST 
does  not  delineate  the;  latitude  circle  at  that  point  but  rather  the  great 
circle  tangent  to  the  latitude  circle.)  The  interior  is  also  taken  to  be 
a, meaningful  section  designated  as  EAST  -f  1,.  Meshes  may ‘take  a 
variety,  of, interesting  forms  from  simple  pie-shaped 'sectors  with- no 
interiors  (posts  all. at  the  same  location)  ,  or  a  representation  of  the 
NORAD  defense  sectors,  or  even  an  inward  spiral  forming  several  layers 
in  an  ingenious  state -like  map.  The  information  about  barriers  and 
meshes  is  converted  and  stored  in  vector  component  form  for  more  rapid 
utilization  by  other  programs,. 

With  proper  key  indications  ATLAS  will  write  the  constructed 
facility  record  on  tape,  one  (as  the  next  logical  record) .  One  may  also 
have  ATLAS  retrieve  a  previous  written  .facility  record  from  tape  one 
with  or,  without  further  card  extensions .  As  many  as  a  thousand  facilities 
in  up  to  thirty  set  s  may  be  placed  in  one:  record.. 


STRIKE 

The  STRIKE  program  can  perform  a  number  of  interesting 
tasks  related  to  targeting.  The  principal  task  is  to  make  a  reasonably 
coordinated  pre-war  allocation  of  various  por  tions  of  a  mixed  stockpile 
of  offensive  warheads  against  selected  sets  of  targets  tinder  various 
restrictions  that  the  attacker  might  wish  to  impose.  The  end  result  of 
this  consideration  is  the  creation  of  a  warhead  record.  This  is-a  iisting 
for  each  of  the  warheads  used  of  their  generated  aimpoints  with  provision 
for  later  insertion  (by  other  programs)  of  origin  of  launch  and  probability 
of  kill  before,  reaching  target.  The  record  iV  segregated1  by  warhead 
types  and  contains,  its  own  keying  information  so  that  itmay.  be  read 
later  if  stored. 

Warheads  may  be  grouped -into  different  types,  because  of 
actual  differences  in  warhead  characteristics  such  as  yield,  CEP,  etc. 
or  because  t h  attacker  would  prefer  to  use  them  in  different  waves' 
(times)  or  for  different  objectives  or  delivered  by  different  carriers. 

Up  to  two  thousand  warheads  split  into  as  many  as  twelve  different  types 
may  be  considered  in  any  one  allocation  and  placed  in  the  same  record. 
The  attack..or  at  least  this  portion  of  it,,  is  described  by  specifying  the 
total  number  of  warheads  of  each  type,,  the  order  in  which  the  itypes  are 
to -be  allocated  (not  necessarily  the  order  of  arrival),  their,  expected 
probability  of  abort  and  CEP. 


STRIKE  presumes  that  there  is  an  appropriate  facility  < 

.  ,  <  - 
record  ih,memory  containing  the  sets  desired  to  be  targeted,.  The 

chosen  (sets  are  request '■'d  iii  any  desired  qrder-.and  for  each  target 

set  -  warhead  type  combination  the  attacker's  desired  constraints 

are  stated.  The’ main  restriction  is  the  most  number  of  warheads 

of  this  type  that  will  be  allowed  to  any  one  facility  of  the  set.  A  blank 

or  zero  here  means  no  interest  in  this  combination  and  the  program 

proceeds  on  to  those  combinations  where  positive  allocations  are 

permitted.  Where  and  within  the  stated  limit  allocation -is.  possible, 

two, basic  choices  are  open  for  determining  the  relative  allocation 

V*  , 

amongst  the  contending  targets  -  requirement-targeting*  or  value 
targeting.  .  . 

If  the  requirement  option  is  chosen,  the  attacker  is  permitted 
to  specify  the  least  number  of  warheads  he  would  like  allocated  to  any 
target  of  the  set.  If  least  is  the  same  as  most>th  \  user  has  completely 
determined  the  allocation  for  this  combination  regardless  of  other 
factors  except  the  stockpile  size,  If  a  tolerance  exists  between  the 
least  and  most  ’boundsj?  third  number  -  a  required  probability  of  kill 
for  each  .facility  in  the  set  -  resolves  the  allocation  to  within  fractional 
-warhead  round  off  (which  may  also  be  controlled)...  The_;pr,Qt 'ability  S>£ 
kill  level  is  however  the  weakest  of  the  constraints  and  will  not  be  held 


to  when  it  lies  outside  the  stated  bounds. 


When  targeting  by  value  options  are  permitted,  these  are  , 

,c  ''  <v  o 

considered  only  after  the  particular  warhead  type  has  fulfilled  its  1 

requirement  obligations.  The  remaining  stockpile  is  then  available 
for  distribution  by  value.  The  object  is  to  maximise  the  total  value 
destroyed.  The  allocations  of  previous  warhead  types  le^.e  the 
target  complex  with  a  reduced  set  of  (expected)  values.  Targets 

■  Jo 

within  the  currently  allowed  sets  compete  with  each  other,  whether 

of  the  same  set  or  not,  based  both  upon  their  remaining  values  and 

their  respective*softness  to  the  current  warhead  type.  .'Softness7  ~  -  ■■ 

is  computed  for  each  potential  target,  as  the  a.verage  single  shot 

probability  of  damage,  based  upon  a  specified  matrix  of  weapon  radii 

for  each  warhead  type  -  facility  set  combination,  and  previous  specifications 

„  o  • 

for  warhead  CEP,  probability  of. abort,  and  target. sizes.  If  desired; for- 
certain  warhead  types,  the  attacker 's/estimaces  of  the  relative  pene-  <  , 
trabilify  of  the  defense  to  specific  targets  can  b  i  introduced  as  an 
additional  modifier  of  the  softness;  The  value  returned  then  is  the 
product  of  the  target's  remaining  value  and  its  softness.  In  maximizing  * 
the  value  returned  for  the  stockpile  available  only  integral  distributions 
of  warheads  are  permitted.  The  mathematics  underlying  the  program  for 
finding  the  allocation  has  been  given  elsewhere*. 


*  Academy  for  Ihterscience  Methodology,  AIM  65-T-I8j  Targeting  Program 
Revised,  by  N.  H.  Painter,  December  1964,  pp,.  5-10.. 


Once  the  allocation  of  the  warheads  is  completed  the  list 
of  aimpoints  Is  constructed.  For  targets  small  in  area  the,  .  ,  : 

aimpoints  coincide  with  the  target  coordinates  and  are  taken  from  the 
facility  record.  For  each  missile  block  the  appropriate  array  is 
generated  from  the  ATLAS  parameters  and  the  aimpoints  spread  out 
over  the  array  in  sequential  fashion  as  far  as* the  allocation  allows. 

For  large  area  targets  ainipoint  patt  'ns  for  each  warhead  type  are 
computed  to  "cover"  the  target  as  uniformly  and  multiply,  as  the 
allocation  permits.  -  < 

For  some  purposes  it  may  be  desirable  to  preserve  the 
warhead  record  as  it  now  exists  in  the  form  of  unaborted  aimpoints. 

This  may  be  done  by  using  blanks  for  the  final  two  STRIKE  data  cards  and 
setting  the  keys  so  STRIKE  will  write,  the  record  on  tape  3.  An 
alternative  would  be  to  proceed  and  read  in  a  new  set  of  probability  of 
aborts  and  CEP's,  which  may  be  the  same  or  different  from  the 
previous  estimated  <ohes.  Jn  this  case  STRIKE  will  convert  the  aim- 
points  to  impact  points  with  the  prescribed  randomization  and  will.  . 
initialize  each  warhead's  probability  of ikill  as  starting  out  with  the 
specified  probability  of -abort.  If  desired  STRIKE  will  immediately 
Morite  Carlo  this  set  of  probabilities  into  rand oih  go-ho-;go  choices. 

The  warhead  record  in  its  completed  form  may  then  be  preserved 
on  tape.  In  either  case  whether  the  record  is  put  on  tape  or  not  it 


remains  in  memory. 


STRIKE  might  algo  be  called  in  order  to  retrieve  one  dr 

more  warhead  records  from  tape  as  needed  thus  obviating  the  need  t 

<  < 

for  reinputting  and  recomputing  attacks  which  can  be  used  multiple 
times.  When  this  kind  of  call  is  made  the-  option  exists  of  either 

r 

accepting  the  record  as  is  or  setting  up  a  new  abort  and/or  impact 
pattern'.  The  latter  would  be  a  particularly  convenient  manner  of 
performing  Monte  Carlo  sensitivity  analyses. 

Finally  one  may  call  STRIKE  again  in  the  same  run  after  a 
certain  amount  of  damage  assessment  has  been  recorded  on  a  facility 
list  and  request  additional  targeting  be  generated.  The  updated  values 
of  the  targets  selected  will  be  used  (whether  the  damage  was  of  the  direct 
or  colocation  nature)  instead  of  the  pre-war  values;  This  presumes  of 
course  that  the  attacker  has  gained,  knowledge  of  this  damage  sufficiently 

O 

well  and  in  time  to  incorporate  it  in  this  targeting  round.  Since  a  hew. 
warhead  record  is  hereby  generated,  the' offense  force  structure  ~ 

involved  may  be  more  of  the  old  force  or  contain  hew  types  or  both. 


Thus  as  one  progresses  from  one  warhead  record  to  the  nextj  references' > 
to  warheads  of  type  one  and  two  for  example  may  retain  their  meaning 


while  types  three  and  four  have  changed. 


-  * 


The  MOTION  program  relates  sets  of  facilities  and  warheads 
to  one  another  in.  various  ways.  In  doing  so  it  makes  extensive  use  of  the 
barriers  and  meshes  which  have  been  set  up  by  ATLAS.  Most  of  the 
global  geometric  relationships  are  explored  by  means  of  the  MOTION 
program.  Besides  generating  information  for  its. intrinsic  value 
certain  specific  information  may  be  useful  in  DEFEND,  particularly 
for  area  defense  systems.  Four  distinct  functions  are  currently  included 
within  the  program,  each  explained  in  the  following  paragraphs.  None 
are  restricted  by  frequency  or  order  of  utilization. 

By  selecting  a  particular  mesh  and  naming  other  facility,  sets, 
the  sites  within  each  named  facility,  are  located  according  to  the  sector 
containing  their  latitude  and  longitude.  A  count  of  the  sector  distribution 
of  sites  is  made  and  printed  along  with  a.  listing  of  the  .sector  number 
for  each  site.  The  seetor  designation  for  >the  site  in  the  facility  record 
is  hereby  tagged  with  the  latest  determination.  Facilities  which  are  not 
found  in  the  facility  record  currently  in  memory  a.re  ignored ;  If  the 
mesh  is  not  in  record,  return  is  made  to  BEAM. 

By  selecting  a  particular  mesh  and  naming  certain  facilities 
which  are  suitable  for  launch  sites  for  selected  warhead  types,  each 
warhead  in  the  sets  chosen  will  be  tagged  in  the  warhead  record  with  an  v 
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appropriate  launch  site.  The  method  for  determining  the  launch  site 


starts  by  locating  the  sector  numbei..  of  the  named  mesh. that  contains 


the  aimpoint  or  impact  .pomt  of  the  warhead.  The  sector  number,  is  then 


1"  i! 


used  to  find  a  corresponding  site  in  the  launch  set.  If,  for  exampl^,  a 


warhead^  fell  in  the  third  sector,  the  third  launch  site  of  the  specified..  , 


set  would1  be  presumed  to  be  the  desired  and  appropriate  site.  There 


should  be  atleast  as  many  sites  in  the  launch  sets  as  there  are  sectors; 


■  >  % 


however,  a  given  site  may  be  listed  more  than  once,  arid  a  given  set 


may  be  used  for  more  than  one  warhead  type.  — Thisjmethod  is  parties 


ularly  attractive  when  range  is  not  the  dominate  routing  criterion  and 


the  user  wishes  to^do  his  own  broad  brush  positioning  of  launch  or 
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penetration  points.  If  the  warheads  are  to  be  defended  against  by  sector 


oriented  defenses,  this  should  be  the  method  for  sector  tagging  the  war¬ 


heads  using  the  same  mesh  as  tx'“  defense.  Otherwise,  the  mesh  need  hot 


correspond  and  could  be  specially  designed  for  the  routing  topology. 


Without  using  a  mesh  but  naming  launch  facilities  for  other 


warhead  types  together  with  maximum  range  constraints  for  each  type., 


a  different  kind  of  selection  of  launch  sites  will  be  made.  In  this  method 


the  launch  site  which  is  closest  to  the  impact  point  from  among  the  de¬ 


signated  facilities  is  found  by  computing  the  great  circle  distances.  If 


the  range  is  withmthc  allowed-  max  this  association  is  accepted,  It  not, 


the  particular  warhead  is  given  an  artificial,  "outrof  -range"  kill  and  will 
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not  be  seen  by  later  defenses  or  cause  any  damage  assessment.  The 
out-o£-range  feature  can  under  certain  circumstances  be  a  planned 
phenomenon  for  achieving  a  modified  targeting  plan.  Short  range  weapons 
used  for  hold  down  or  defense  roll  up  might  well  be  targeted  this  way. 
Counts  of  the  out-of-range  kills  as  well  as  the  loads  placed  upon  the 
penetration  or  launch  points  in  order  to  meet  the  in  range  targets  are 
printed  out. 

By  naming  appropriate  barriers  for  selected  warhead  types, 
for  whom  launch  sites  have  been  previously  designated,  information 
relative  to  the  attack  crossing  these  barriers^can  be  obtained.  For 
each  warhead  the  position  where  the  great  circle  arc  from  lairnch  to 
target  crosses  the  selected  barrier  is  found.  Not  every  warhead  may 
cross.  For  some  the  target  may  lie  outside  the  barrier.,  i.  e. Honolulu, 
while  others  may  be  launched  interior  to  the  barrier.  Other  paths  may 
cross  more  than  once  in  which  case  the  first  computed' crossing  (lowest 
numbered' arc)  is  used.  For  each  arc  of  the  barrier  the  sum  of  the 
probabilities'  of  survival  of  the  penetrators  crossing  that  arc  (surviving 
penetrators  if  Monte  Carlo  has  been  made)  is  printed  out,  as  well  as  the 
weighted  average- of  time  to  go  from  arc  to  target.  This  information  is 
not  however  currently  retained  in  memory. 
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The  DEFEND  program  is  called  for  the  single  minded  purpose .- 
'of  increasing  the  attrition  on  the  attack.  The  program,  may  be  called  at  -  f 
various  times  and  in  each-call  selected  sets  of  defense  facilities  may  be 

■  •  .  ,  C,__ 

requested  to  do  battle  in  arbitrary,  order  or  repetition.  A  variety  of  local  • .  •• 
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area,  missile,  and  aircraft  defenses,  may  be  included.  In  each  defense  ■  A,  n  .  „  o0  vj 
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battle  tire  defense  facility  is  faced  with  the  presence  of  a  possibly  mixed  „  a  ...  i>  '>  s 


selection  of  offense  types  in  a  common-time  span.. 

Several  geographical  mechanisms  are  available  for  confining.  ' 

c  O 

v\  the  general  commitment  Of  each  defense  site  or  sector  .  Within  the  0 
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confined  environment  the  efficacy  of  the  site's  normal  attrition  capability 

< 

becomesAimitcd  with. increased  rapidity  of  events  as  the  total  raid  is  -  . 

compressed  in  time.  The  actual, locale  of  the  battles  are  not  sought  nor 

*  C 

are  detailed  intercept  assignments  or  geometry  computed.  Given  reasonable 
computation  of  forces  which  the  defense  can  mount  and  the  segments  of  the  < 
offense  they  are  directed  at  it  is  not  unreasonable  to  manage  the  intercept 

•  c  * 

variations  as  a  random  application  of  kill  potential  to  moving  targets.  » 

„  ,  * 

Each  defense  site  listed  in  the  atlas  facility  takes  its  tvirnat  < 

doing  battle  in  order.  The  simplest  commitment  mechanism,  is  tha.t  appro -■ 
priate  to  local  terminal  defenses.  The  sites  have  a  circular  radius  of 

i 

coverage  around  their  physical  coordinates.  This  coverage  may  vary  ,  •  .«* 
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with  penetrator  types  including  'Shrinking  to  nothing  at  all.  Those  warheads 
which  are  scheduled  to  impact  within  the  coverage  and  which  are  still  alive 
at  the.  time  the  site  is  ready  to  commit  constitute  its  threat.  Note  the 
coordination  which  is  thereby  achieved  by  overlapping  sites  since  a 
first  site  kill  will  not  be  a  second  site  target  but  a  first  site  failure  Will. 

A  second  commitment  mechanism  is  more  appropriate  to  area 
missile  defense  as  such  sites  are  not  likely  to  be  deployed  near  targets 
they  defend.  In  this  an  effective  circular  intercept  coverage  around 
the  .sites  are  assumed  as  before  but  this  time  at  some  altitude.  The 
footprint  or  ground  projection  of  the  equivalent  coverage  .depends  upon 

"•  '  or 

the  direction. and  angle  of,  reentry  of  the  missile  attack.  With,  such 
information  for  the  various  missile  types  the  down  range  displacements 
are- inserted.  A  similar  type  coordination  occurs  when  footprints  overlaps 

-  v 
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Either  of  these  geometrical  partitioning  may  of  course  be 

-  •  0 

applied  to  aircraft  defenses,  particularly  of  the  Hercules  or,  Bomarc  type. 
For  long  range  fighters  the  facility  ‘bites M might  be  taken  as  centers  of 
NORAD  sectors  of  responsibility  or  air  cap  positions  displaced  up  range  to 

v 

the  center. >of  battle  coverage,,  although  the  “reentry,"  is  awkward  in  the,  latter 

*  w  o 
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case.  A  third  .procedure  is  to  request  defense  by  sector.  This  entails 
that  a  mesh  appropriate  to  the  sector izing.be  included  in.  ATLAS,  the 
defense  facilities  located  by  sector  by  MOTION,  and  the  defendable 
:penetrator?  sector  tagged  by  MOTION  with 'the  same  mesh.  Each, site 
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then  battles  those  penetrators  with  the  corresponding  sector  designation. 

A  fourth  procedure  is  to  use  a  simultaneous  constraint  of  radius  of 

■* 

coverage  and  sector  coincidence,  but-an  apjAica-tion  is-not  apparent. 

Whatever  the  commitment  partitioning  procedure  is,  an 
actual  count  of  surviving  penetrators  is  made.  For  the  defense  environment 
point  of  view  this  count  is  modified.  If,  for  the  span  of  time  allocated  to 
this  raid,  a  group  of  penetrators:  can  be  presumed  to  be  all  in, the  battle^ 
area  and  remain  there  throughout  each  then  has  the  same  dwell  time 
as  the  raid  time.  If' a  second  group*  of  fast  moving  penetrators  all  pass 
through  the  battle  area  during  the  raid  interval  it  is  assumed  that  they  enter 
ahu  leave  at  random' tunes.  Each  has  a  dwell  time  that  is  only  a  fraction 
of  the  raid  time  and  the  number  present  at  any  moment  will  be  on  the, 
average  this  same  fraction.  A  third  group  of  penetrators  may  enter 

and  leave  the  battle  area  at  random  moments,  not  because  of  their  speed,  ; 
but  because  their  entries  are  straggled  out  over  a  long  time  relative  to  the. 
raid' interval.  Their  dwell  time  should  also  be  shortened  to- reflect  the 
lower  density  present  during  the  raid.  However 4  when  handled  this  way 
the  defense  only  receives  justice  if  repeated  raids  are  called  on  the  stragglers. 
Generally  little  is  lost  computationally-/  by  time  collapsing  a  straggled  raid, 
into  one  separate  simultaneous  raid  if  -saturation-and  inventory  ready 
allowances  are  made.  .  .  v- 

■v 

The  count  of  penetrators  present  may  be.  further  reduced  by  a 
visibility  factor  or  increased  by  u\e  portion-of  the  cloud  of  false  bogies 
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accompanying  the  visible  penetrators  which, cannot  be  successfully 
discriminated  against  by  the  particular  defense;  This  total  density 
of  bogies  persist  throughout  the,  raid  time  a^hough  for  fast  penetrators, 
the  identities  may  be  changing.  The  hiimber  of  tracks  a  site  or  sector 
is  able  to  maintain  generally  has  an  upper  limit  .so  that  the  probability 
a  given  bogy  is  a  track,decreases  as  the  number  of  bogies  approach  the 
saturation  limit.  This  probability  of  track  is  taken  of  the  form 

-1/  BOGIES  \2  • 

PT  =  e  2  ^SATURATION./ 

If.  the  bogies  greatly  exceed  the;  limit  the  number  of  actual  tracks  may  in 
fact  decrease;  Each  interceptmade  oh. a  track  requires  a  ground  data 

*  v 

link  or  communication. channel,  an  interceptor  salvo  and-a  certain 
amount  of  time.  The  total  number  of  intercepts  made  during  the  raid 

t  *  • 

may  be  limited  by  channel  duty  rates,  interceptor  ready  inventroies  or 

J 

failing  either,  of  these  a  command  limitation  when -a  satisfactory  number 
of  intercepts  per  distinct  bogy  is  attained.  (If  the  average  dwell  time  is 
half  the  raid  time,  the  number,  of  distinct  bogies  will  be  twice  the  number 
thatare  on ‘the  average  present.)  The, interceptor  inventory1 may  or  may 
not  be'  i*equesJ  d  to  be  counted  down  as.  intercepts  aremades,  There  are 
■  r'aflphal.appiicatibhs  for  both  missile  and  fighter  type  defehses.ior  either 
request.  The  number  of  intercepts  per  track  per  unit  time  is  computed. 
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The  probability- of  survival  of  every  warhead -within  the  coverage.decreeses  ^ . .. 
in  relation  to  the  probability  of  kill  of  the  pehetrator  type  per  intercept^  ;V  ‘  '  ^ 
the  number  of  intercepts  per  dwell  time,  and  the  probability  of  being  tracked.  -  ^  5 

*  »  ~~  4'  „ 

If  desired,  the  expected  probabilities  of  kill  may  be  left  to  accumulate  ,  -=^r 
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through  several  battles  or- they  may,  be  Monte  Carloed  each  time.  It-ds',;  ' 

generally  desirable  to  ha  ve  each  warhead- Monte  Carload -at  least  the  • 

c% 

last  time  before  entering  the-  damage  assessment  program.  The  expected 
number  of  defense  kills  and  opportunities  for  each  site  and  penetrator  type 
is  printed  out  for  each  battle.  ’  * 
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The  ASSESS:  program  manages  all  .the  damage  assessments 
desired  on  any  or  all  real  facilities  whenever  warheads  have  passed  '  -  . 
through- tho  dafohsoi  Facilities  dp  hot  have  to  have  been,  targets  and  ■ 
colocation  damage  is  accounted  fori.  When  it  is  called  into  memory, 
by  the  'BEAM .-monitor,  the  next  four  input  cards  are  read.  These  cards 
•Establish  the  desired  yields  (in  megatons),  the  thermal  efficiency  i 
adjustments  (normally  unity) the- fission -fusion  fractions,  and  the 
altitude  type  of  burst  (surface,  air  or  high  altitude)  for  all  of  the 

«  i 

various  warhe? 4  types.  The  program  contains  numerical  data  which 
define  blast  contours,  thermal  contours,  and  fallout  contours 7  This 
basic  data  is  for  a  10'MT  warhead  with. 50-50  fission-fusion.  For  bias* 
these  contours  are  concentric  rings  of  lbs  per  square  inch  overpressure 
at  various  nautical  mile  distances  from  ground  zero  (Figure  1)  For  „ 
thermal  they  are  rings-  of  incident' calories  per  square  cm  at  various  ■ 

v 

nautical  mile  distances- (Figure  2).  .For  fallout  the  contours  .are- a  ‘ 
nested  set  of  oblong  rectangles  .skewed  toward  the  down -wind  (easterly) 
direction,  for' total  roentgens  received' in  48  hours  with.no  protection 
(Figure  3).  The  blast  and  thermal  contours. .change  for  each'  of  the 
three  di-ffercni-altituue  type  bursts,  while  the  early  fallout  is  assumed 
to  be  appreciable  only  for  surface  bursts-.  For  a  particular  warhead 
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type  the  appropriate  altitude  contours  are  scaled  to  .the  correet,  .  *<?  v  *• 
yield;  for  blast. and  fallout  the  scaling  is  by  yield  to  the  1/3  powers  ,  .• 

for  thermal  by  yield;  to  the  .-42  power.  ■  ’ 

•  '  ^  <■  .r 

A  particular  set  of  warheads  arc  withheld  from  considera.ti.on 
during  the  current  assessment  if  their  altitude  types  are  set  to  a 
negative  one..  This  permits  different  offensive  components  of  the  same 
strike  to  be  assessed  by  different  assessment  calls.  One  advantage,  of 

r 

this  procedure  occurs,  when  the  strike  contains  some  time  delayed  r ; 
elements  and  the  defense  survival  and  effectiveness  can  be  studied -.and 
altered  between  waves.  An  alternate  procedure  which  is  some  times,  t 

A  V  c,'  •  2> 

useful  is  to  construct  one  "strike  .kit- 1  such. that  d iff er ent  c ombinatipns 

of  warhead  types  represent  different  threat  levels  or  targeting  objectives. 

Each  assessment  then  represent  the  effectiveness  of  one  of  the  threats. 

In  one  application,  the  low  threat  was  obtained  by  taking  only  warhead 

•  * 

types  1  -  5,  the  medium  threat  was  obtained  by  augmenting  the  low 
threat  with  those  warheads  called  type  6,  and  the  high  threat  by  the 
further  augmentation  of  warhead  type  7.  Warheads  of  type  5,  6,  and  7 
were  all  physically  alike,  the  distinction  being  only  in  the  extent  of 

f 

targets  hit.. 

*1 

When  the  warhead  initializations  are  complete, 'requests  to 

<r  v 

assess  damage  to  classes  of  facilities  are  processed.  Each  subsequent/  . 
input  card  requests  one  class  in  the  facility  atlas  be  assessed  and  „ 
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supplies  the  levels  of  blast,  thermal,  .and  fallout  stress  (or 
protection)  considered  to  be  acceptable  for  that  class.  Oh  the  basis 
of  the  class  name:  supplied  (and  .previous  ATLAS  information),  the 
subroutine  distinguishes  the  heed  for  five  different  .procedures, 
depending  on  whether  the  class  represents  1)  bases  or  sites  with 
similarity  in  area  and  function,  2)  strategic  missile  squadrons  or 
"blocks”,  4)  urban  areas  with  shelter  data  (on  tape  prepared  by 
the  SHELL  subroutine),  or  5)  "rural  population'1  or  continental  map 
of  values.  Most  military  installations,  with,  the  exception. of 

i 

missile  blocks,  are  usually  considered  to  be  treatable  by  the  iiri  V  ‘ 
procedure,  and  this  is  the.  one  discussed  in  detail.  The  area  of  the 
base  (zrR^)  times  the  density  of  evaluation  points  per  square  nautical 
mile  (both  inputs)  determine  how  many  evaluation  ,points  are  to  be  - 
placed  around  each  base.  This  number  is  at  least  one  but  no- more 
than  400.  Whatever  the  number,  a  squarish  pattern  of  that  many 
points  is  constructed,  centered  on  each  base  and  scaled  to  cover  the 
equivalent  area  of  the  base  (Figure  4). 

For  each  base  in  the  facility  class,  the'  list  of  warhead 
impact  points  is  examined.  Warheads  which  have  aborted',  or  been 


successfully -defended  against/  Or  'disallowed 


mentioned  above)  are  skipped.  Those  impacts,  which  are  too  remote 
to  affect  the  particular  bases;  are  also  skipped.  For  each  remaining 
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warhead  the  appropriate  effects  contour  s  are  placed  around  the  ,  -  . 

impact  point.-  The  location  of  each  of  the  base  evaluation  points  is: 
found  with  respect  to  the  three  sets  of  contours  and  the  respective,  .  ,  o 
PSI,  calorie,  and  roentgen  doses  that  would  be,  received  from  each 
detonation  at  the  evaluation  point  are.  weighed.  For  blast  and; 
thermal  doses  the  point  is  charged  with  receiving  .the  maximum  levels 

«. 

x  V 

found  among  all  the  warheads  detonating  while,  for  fallout, the  dose  is  ,  0' 

the  cumulative  of  each  warhead's  contribution.  In  order  to  present  0; 
<to  the  analyst  an  indication  of  the  effects  experienced  on  each  fac ility  .  ^ 
of  this  type,,  the  three  types  of  dosages  are  averaged  oyer  all  the  . 

■> 

target's  evaluation  points  and  these  averages  are  pointed  for' bias ti 

„  >\ 

thermal  and  fallout  before  target  vulnerability  to  the  "absorbed"  effects, 

are  considered..  :  " 

”  < 

When  the  doses  have  been  computed  for  all  the  evaluation 
points  around,  the  base,  the  base  survival,  is  considered..  The  probability  , 
a  particular  evaluation  point,  considered  to  represent  a  portion  of  the 
base  ,  survives  is  determined  from-  functions  of  .the-  effective  .doses  - 

wi  - 

received,,  which  are;  the- actual  doses  relative  to  the  {iupsit)1  stress  - 

.  v,  - 

limits  for  the  facility  class.  For  blast  and  thermal  the,-  functional 

2  ..... 

relation  is  e  ‘{Figure.  5),  while  for  radiation  it  follows'  a  fatality  ■ 
table  (Figure  6) .  The  fraction  of  the  base  as  a  whole  which  is  killed 
by  blast  is  the  number  of  evaluation  points  \yhich  do  hot  survive  the  -  ,V 
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the  blast  alone.  The  fraction  of  the  base  killed  by  thermal  is  the  ? 

*  1  ' 

count  of  points  which  survived  the  blast  wave  but -did  riot  survive  the  •'* 
fire  or  thermal  effepts.  The  portion  of  kill  attributed  to  radiation  ‘  ■>  . 

would,  ba  those  which  §§e&ped  the  ftees  two  etifect©  y9t  received  a 
fatal  roentgen  dose  in  spite  of  the  protection  factor;  To  find  what  ■' '  ■- 
effectiveness  fallout  alone  would  have  an  assessment-call  could 

,i„  "V 

be  iriade  where  the.  other  effects  are  "turned  off"*  The  value,  of  the 
base  is  assumed  to  be  uniformly  distributed  over  its  area  or  equivalently 
its  points  and  tho  value,  surviving  that  is  printed  out  as  that  portion  which 
survives  all  three  effects.  If,  as  sometimes  occur.,  multiple  independent 
assessments  are  desired,  the  option  of  leaving  the  base  value  in  ,  . 

memory  unchanged  may  be  specified.  Normally,  the  value  decreases 

V  v 

with  succe.  sive  assessments.  The  total  values  surviving  or  killed  by 
each  of  the  effects  is  summed  and -printed  out  for  each  facility  class. at 
the  end  of  its  processing,  A  blank  card  ends  each  assessment  call-arid 
returns  control  to  the  BEAM  monitor.  ‘  ; 

The  rural  population  does  not  occur  as  a  facility  in! ATLAS.. 

It.  or  any -other  regular  geographical  array  of  values,  such  as  rural 
MVA,  may  be  assessed  by  directly  reading  iri  a  density  "map";  Each 
1°  by  1°  latitude -longitude  square  corresponds  to  a  base  prurbania.rea- 
in  the  normal  procedure  .and  has  its  own  value.  A  specified  square  . 
density  of  evaluation  points  are  arranged  symmetrically  to  cover 
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every  cell  (max  10  by  10)  and  the  nation  has  a  regular  pattern  of 

__  $ 

evaluation  points  except  for  the  adjustment  which  is .$&.&»  -for  the 
cosine  of  the  latitude.  The  value  of  a  cell  is  represented  as  an  alpha-, 
numeric  character  oh  a  scale  0-50  with  an  arbitrary  scale  factor  ior 
the  map  as  a  whole.  Each  input'  card  gives  the  latitude  string  of, 
characters  for  a  given  latitude  band,  and  the  set  of  cards  are.to 
latitude  order  north  to  south.  Blanks  cells  are  not  assessed.  The 
ouiput  is  represented  in  the  seme  format  as  the  input  with  maps  for 
the  survivors,  blast  km.  and  fallout  kill  distribution.  Cumulative  '• 

nationwide- totals  are -also  printed  out. 
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SHELL. 


This  program  is  concerned  with  providing  data  on  the  .  '' 

blast  and. fallout  protection  'the  population  might  have  in  various  <  l 
circumstances.  The  output  can  be  made  available  to  ASSESS  which 

< i  1  ’• 

enables  the  latter  to  perform  a  more  sophisticated  evaluation  than  o 
it  would  otherwise.  Three  distinct  functions  can  be  performed  using 
SHELL,  '  ' 

A  basic  set  of  data  is  required  to  describe  the- physical 

-•  '  <r 

number  and  nature  of  the  shelter  spaces  that  exist  at  each  facility  site 
The  natural  facility  is  of  course  the  urban  areas  and  this  is  the  one  -  . 

for  which  a  data  bank  has  been  accumulated  and  used.  There  is  no  ., 

•> 

reason  why  similar  data  pould.not  be  applied  to  military  or  other  ..  ^ 
facilities  that  contain  special  groups  of  people.  The  national  shelter 
survey  was  the  source  for  creating  the  urban  shell  deck.  Sinc.e  the  * 
origin  form  of  the  data  was  not  accumulated  in  a  manner  to  be 
applicable  to  urban  areas,  selected  component  data  was  punched  on 

O 

cards  for  various  parts  of  the  urban  complex.  This  consists  of  the 
number  of  spaces  in  the  standard  PF  and  PY  categories,  defined  by. 
OCD.  The  first  function  and  use  of  SHELL  is  to  tabulate  and  order 

,  -  -  v  *  '  '  ’  _ _ „ _ c 

the  totals  to  match  the  site  list  in  the  facility  record,  in  this  base 


the  213  U.  S.  urban  areas  arranged  by  rank  (identical  to  the 


Certsus: 


Bureau  ranking).  This  is  one  instance  where  order  in  the  facility  record 


is  important.  The  tabulation  is  usually  made  once  since  survey  data0 
generally  do  not  change  often  and  a  permanent  shell  tape  written.  .  , 
Thereafter,,  the  master  tape  is  used  to  generate  the  Operating  tapes  / 


{or  SHELL  and  ASSESS.  Both  assume  that  the  urban  deck  matches' 


in  number  and  order  of  places.  , 

The  second  function  that  SHELL  may  do  is.  to  construct 
and  add  additional  shelter  spaces  as  directed.  The  new  spaces  may  be 
constructed  at.a  specified  PF  or  PV  category  or  both;.  The  number 

’  .1 

-  •  e 

added  to  any  given  place  will  depend  on  if  and  how  much  deficient  that 

place  falls  below  a  requested  criterion.  The  implication  in  terms  of 

,  > 

« 

the  national  total  spaces,  constructed  is  printed  out.. 

The  third  function  is  to  distribute  at  each  place  fts  popu» 


lation  value  in  some  manner  into  the  various  shelter  spaces  it  has 
available..  One  choice,  which  might  be  suitable  when  long  strategic 
warning  occurs,  would  be  to  take  maximum  protection.  In  this  case ' 
the  hardest  spaces  are  filled  until  either  there  are  no  more  people  or 
no  more  spaces.  An  opposite  or  no  warning  case  would. leave  a  .  • 
percentage  of  the.  population  with  mmimum  protection;  The  remainder; 
fill  up  the  spaces  in  a  pessimistic  fashion,  first-next  poorest  and  on 
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up  perhaps  leaving  the  best  spaces  empty.  The  most  realistic  1 
option  assumes  that  all  shelters  regardless' of  type  are  filled  to  a 
certain  percentage  of  capacity,  the  percentage  utilization  specified. 
Any  remaining  population  has  minimum  protection. 

Any  number  or  mixture  of  utilization  can  be  stacked  on 
a.  temporary  6r  permanent  tape  for  concurrent  or  later  ASSESS 

*  *  v 

requests.  ‘ 
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SURVEY  . 

The  SURVEY  program  has  been  brought  under  the  operating 
command  of  the  BEACON  monitor  and  ib  exercisable  *  II;  is  currently  ■ 
more  autonomous  \d  thereby  less  in  the  battle  stream  than  the.other  , 
programs.  Wh<*  ie  Task  41  study  settles  upon  an  MCI. attrition 
model,  it  potentially,  may  relate  strongly  to  SURVEY  and  together 
integrate  more  firmly.  At  present  SURVEY  may  be  tisecl  to  explore 
questions  concerning  missile  and  satellite  trajectories,  their  visibility 
by  one  or  a  group  of  sensors,  and  potential  warning  or  intercept  aspects 

A  number  of  objects  may  be  launched  from  arbitrary,  places 
and  directions  around  the  globe  with' simple  or  elaborate  descriptions 
of  the  boost  stages' and  thrust  forces.  The  gravitation  attraction  forces' 
by  the  oblate  earth  are  accounted  for  and  the  drag  forces  change  with 
altitude  and  cross  section  in  a  density-temperature  varying  atmosphere . 
The  spatial  and  rotating  earth  trajectories  are  integrated  from  the 
force  equations.  A  trade-off  between  speed  (of  computation)  and 
accuracy  may  be  made.  •  • 

A  group  of  sensors  may  he  deployed  any.  wherev  These 
sensors  which- contain  the  trajectory  in  their  horizon  are  indicated’ 
as'  a  function  of,  time  as  well  as.  the  object  position  in  sensor  coordinates 
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COM  AND 


;>c  ..  ^  H-  -  .* 

'  *'  ^ 

,<  A  »  '*  -  *  <  A% 

-i 

"  -•  ,:'  ;;i 

*'  *  ^  4  Vi 


c  v  *0}  t  *  v,  *  ' 

T.he  GOMAND  program- is,  a  .utility  aid  in-  providing  visual?,/,*  ■«- '  '  : 

-  ~  *  -  *  &  0 '  ’  „ 

Jr^/i  .  * 

plots  of  various  geographical  relationships  between-  .either  sets  in  the;/  -\*y 

*  *  /T. 

facility  records  or  in  the  warhead  records, or  both  simultaneously.  *  * 

A  separate  map  is  produced  with  each  call.  With  a  1  °  ' c e,ll  r e s oluti on  ,  / 

,  '*  <>  ^  h’  <> 

the  map  frame  can  be  as  large  as*  50°  in  latitude  by  120°  in.  longitude*  -  v-i**:  ■  - 

■  ,  '  -  Ci  __  __ 

Larger  or  smaller  areas  of  the  globe  may  be  framed  with  a  corresponding- 

*  * 

loss  or  gain  in  resolution.  The  frame  may  be  positioned  anywhere',*  The  i\ 

.  o 


A 


U.S.  may  be  nicely  framed  by  a  25°  x  60°  map  with  Vji  .cell  ..sizes,-  '• 


\  >  *.s  , 


By -naming  a  record,  facility  or  warhead  list,  a  facility  set^ 
or  type,  and  supplying  a  symbol,  the  program  will  print  that  symbol 
in  every  cell  containg  such  a  site  or  impact.  Several  requests-may. 

.  '*  C.. 

be  made  for  the  same  map  when  pne  wishes  to  yiew  the  correlation 
between  them.  Clutter  however  should  be  avoided. as  multiple-  entries 
to  the  same  cell  result  in  converting:  the  symbols  to  asterisks-./ 
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UTILIZATION  OF  THE  PROGRAMS 


The  BEACON  Monitor  is  initially  given. -program  control..;  ,. 

The  .first  card  of  . the  data  deck  is  a  date  card  which  supplies  the . date*  « 

i  ‘  .  r  *-v  " 

a  name  or  number.,  and  eight  words  for  identification  of  the  run.  This.-  s 
information  is  printed  out  at  the  very  beginning  as  part  of  the  heading 
for  the  run..  The  monitor  is  able  to  read  cards  in  any  format;  it 
passes  them  by  until  it  comes  to  a  BEAM  card,  which  ds.  identified  by 

the  Word  "BEAM*  in  the  first  six  columns  of  the  card  (the  quote  sign  " 

* 

signifying  a  blank  character).  The  next  word  (col  7-12)  is  a  subroutine 
name  and  must  be  spelled  precisely  as  one  of  the  available  routines 
ATLAS",  STRIKE,  DEFEND,  ASSESS,  SHELL",  MOTION;,  SURVEY,  . 
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COMAND  or  be  completely  blank;  Any  other  combination  of  characters, -  ■  ;  , 

such  as  END,  acts  as  an  end  of  run  signal  to  the  monitor.  After  the  -  a  O; 

'  '  '  '  “  '  . '  Z  '  i  ) 

u  <V 

monitor  has  processed  the  remaining  information  on  the  BEAM. card, 

.  c 

control  is  then  passed  to  the  appropriate  routine  requested..  As  each  -•  .  ;  >>  * 

4  <>  -C 

routine  completes  the  work  for  which  it  was  called,  control  is  returned  1  <  *  , 

to  the  monitor  which  resumes  scanning  the.  data  deck  for  more  BEAM,*  . 

* 

cards.  It  is  in  this  manner  that  the  user  achieves  His  participation  ir*  -  -  s> 
the  battle.  By  dividing  the  problem  .into. appropriate -segments,  ^  '  \  .  : 

arranging  their  order  and  inserting. BEAM*  cards  as  definition  of  this  .  " 

.  ,  -  ’  *4 ,  .  ‘  , 

control,  the  user  literally  commands  the  computer  to  perform  the  '.'.V,  ! 

*  =  v  ,  J./  <, 

detailed  calculations  for  a  battle  which  he  outlines.  .  v 
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Examples  of  some  hypothetical  deck  set  ups-iix^y  be 
helpful.  Simple  preliminary  jobs  are  given  first* 

■-  y-!/  ,'1-  - 

•  .  i 

1}  Form  a  raw  shelter  tape  for  permanent  use  from  the 

•  ’  '  i  ♦  "  1  '*  '  C 


survey  data,  no  construction,  no  utilization." 


Date  card 
BEAM*  SHELL1’ 

[PF  and  PV  decks]' 


;2?2 
keys|  5=0 
>  6=0 


Blank  card 
BEAM*  END 


■2)  Run  a  missile  trajectory bast.BMEWS. 
Date  card 
BEAM*  SURVEY 

t  , 

[Sensor,  lisf] 


3) 


.[Launch  cobrd  data] 

Blank  card  " 

BEAM*  END  *  • 

Write  the  main  part  of  an  atlas  on  tape. 


Date  card' 
BEAM*  ATLAS 


iISAG:  base  si 

H  ]•  " 

ft' 


[Nike  X  sited] 

V  ,  -  '  ~ 

■1  -•  38 

'  .  *  - 


Key  2°2  (write  tape) 


:r  bases] 

•'[Sage  sites] 

>  , 

[AWACS  Barriei], 
[NORAD  Mesh] 


\  -  9'-  ..  i 


o'  ‘  ,  '< 

'  6 


i  A-  ^  » 
b  ' 


.  • 


'  3f-  , . 
0-  v  ■ 


w 


E  ank  card 
BEAM*  END 


Read  previous  atlas  from  tape,  add  a  set  of  launch  sites, 
target  and  put  aimpoints-  on.  tape.,  .find  dwelt  times  for 
fast  and  slow  bombers  from  AWACS  barrier  to  target^, per 
NORAD  sectors.  - , 


O  v, 


■ft"  o-  i 


<1. 


£  :> 


Date  card 
BEAM*  ATLAS 

[Launch  sites]- 
Blank  card 
BEAM*  ‘STRIKE 

[Warhead  characteristics] 
:get  requests]; 


Key  2-1  (read  tape) 


f  'V 


* 

ft 


3 


,1  * 


Key  2  -  2 


[estimated  PABORT, 


_<  j. 


2  blank  cards  (aimpoints  only,  ,no  actual  aboits) 


;o 


1-3.9 


<>*«  -  -  ?l 

t 

■’v  _  : 

*5  C  J7 


v  9 


C(  a,,  ✓  ,  *■  «■  b  i*  v  1  -1  -  $C 

b  ^  ^  \  *  <i'  '  >  *  y 


BEAM*  MOTION 


.  f 

f  >  • 


NORAD3 


NORAD2 


NORAD0 


(fast  type)  {'low  type) 

LAUNCH.  LAUNCH  ■{ 


AWACS 


;;  o~/  <.  <*  * 


AWACS'  {fihd.c'^ol’IIhgs/ 
'and^timfesjA«'  * 


Blank  card 


BEAM*  END 


Run  three  stage  defense  battle  and  assessments. 


Bate  card 


BEAM*  ATLAS 


} "  ( 


Key  2  -  1  (read  atlas  tape) 


Blank  card 


BEAM*  STRIKE 


Key  2=1  (read  warhead  tape) 


;; 


BABT  and  CEP 


Monte  Carlo  impacts 


BEAM*  DEFEND 


[Missile  attack  only] 
Nike  X  -  .parameters1 


Blank  card- 


BEAM*  ASSESS 


fix**, _ .. 

I'miooiikv  •ui»wy>vtv  rw***.jr 


SAGE  -  hardness  data. 


BOMARC  -  hardness  data 


■  V  *  •>/?.'■:  .  i  .  t  - 


r  4‘, 


< .  <**  • 


,  *  Qv .  ^  <<y  ?•"  i 


MAP  OF  AREA  BOUNDED  BY  25.  AND  49,,  DEGREES  ‘LATITUDE 


D  • 


'  D 


R’ 


1  D 
DD1 


R 


R* 


•R 


« 

D 


G  R 

«  *• 


R 


THIS  SECTION  IS  CpMAND 

! 

I 

MAP  US  49v&0  >126.00  24.00 


*o  <> 
■  < 


6  • 

,  c 


n  'f<7'  '-VV^ 


£?* 


:  -  .1  ,  ,*  • 

o  ■’ 

<■  <■ 


r  ^  f*-t 

.  >  !'*■, 


.  f  . 


- '  ,«  *-  : 


.  <, 
■'V  , 


c 


R 


D* 


D  » 


RD 

f'» 

SG 

:fg 

UU 

j-u 

RB 

FR 

OANURB 

FC 

59.00  0.30  0,30 


SECTION  II 


t  ^ 

**  '*  t)  c''  -; 


*  *'  O  .*1  , 


,«*  .*<>  ‘  C* 


•\ . .  •?>  •* 


BEACON  MAIN  ROUTINES  for  the  IBM  70,94 


1.  BEAM 


2.  ATLAS' 

3.  STRIKE 

4.  MOTION. 

5.  DEFEND. 

6.  ASSESS. 

7.  SHELL 

8.  SURVEY 

9.  COMANB 


o  ’  >-.t  ’<.* 


c 


oj'v, 


STB  ETC  BEA*\)K  BECK  _ -  -  -  - 


CARDS  rQR  lEArO  MONITOR  PROGRAM 

COMMON  /M3N/  IDENTC8)  4KEY(12)  4LBLC4)  i'NUM  . 


COF'*0‘  p-AR/  *. A « L R » RADI (JS * P I E.VR  ANDA rRANDB 4  '  '' 

t  COMV«SCALAT«S.GALONaPOLAT,POL6N, 


+  ROUn'Di EPSI  » DILATE *CPKl-dCPK2*CPK3 
.  ..  DIMENSION  NEXT  (9)  ,  .'APUN  (  A)  . . .  .  . 


DAT  A  ICHECK/6H  PEAM0/4KEYHOL/6H] OOOOO/4 
.  +  TAPUNM) /24H01,  ,A3+++**+n3.;A4G4;l04/4  _ 


XTfc  1W54HATLAS  STR  IKEDEFENDASSESSSHELt.  .MCTJONSURVEYCOMAND 
NAMELIST  /MIJPAR/  RAN0A-«RAnD8 »LA<»LQ 4  _ ___ 


*  CONV^SCALATiSCALON^PQUAT^POLON,  . ■ 

t  ROUND >  FP Si « D I»L  ATE* CPK>i  o CP.K&VCPK3 «  . .  . . 


COMMENCE  '  * . 

REAn  (LA»A10)  DATE«  JOB 9  ?DEN? . .  ...  - 


WRITF  (LBi400>  IDENT 9  DATE*  JOB  .  *  " 

DO  1  07  L  .  ~  1 4  4  .  ...  .  ...  _ _ 


'  !F  <  L  4  NF ,2  )  REWIND  1  .  “  '  '  ‘  .  ” 

107  „  LPLM  )  a,  0  . . 

101  -PE AD  (LA 4 410)  MATE 4  NAME 4  IDENT«.  KEY 

ir  (  ICHECK  .NE;  MATE V  .60  To  101  .  ..  . 


r  continue  . 

STOP  .  . 


I  DO  104  -K  =  I  4.I2 .  • 

KEY (K)  =  KEY(K)/KEYWOL 


F (  KEY  <K ) WLT . 0  >  KEY  <K)  «  0 
K  CONTINUE  .  ... 


NMM  =  1000^KE.Y(9T  «  100OKEYMQ) 
T  F  (  KEY (2  ) • L  T . 7  )  GO  To  105 


READ  (LA4NUPAR) 

...WRITE  .(LBvNMPAR)  . 


>  lORKEY(ll)  *  KEY C12') 


106:  WRITE 
„  DO  9fl 


(LB  4 40 3  ) 
L  =  1*4 


NAME',  1 0£NT« 


*FrKYE<HO  r  GO'  TO  90 . •  ' 

GOTO  (.  91 4  9  2 , 9 3.»  9  4  0  9  5 « 9  A ,  97  *  9 8  *99)  #  K 

91 

PRINT  414  IDENT(L),TAPUN(II)»!DENT(2) 

IBUU  a  0  .  .  •  .  '  ' 

'  LAPSE  -  L . 

GO  TO  90  .  . 

• 

»  t 

■ 

s 

'  92 

READ  (L)  LBL“(L) . -  '  '  ' 

go  to  90  ...  :  . 

93 

BACKSPACE  L 

LPLd  >  =  MAXo  (L8L  (L )  <-l  «'0) . 

95. 

GO  TO  90 

REWIND  L  . . . . . . .  .  ... 

- Q  1  »  1  t  ...  fk  —  '  '  '•  “  ’  .  '  -  '  ‘ 

.  V  -  * 

TFT  KEY (8) ,E0«2  ) 
IF(N)  84490,82  . 


82  no  83  r  a  r«N. 

83  RE  ATI  CD  LBL(L) 


GO  TO  90 


2-2  ' 


in  LAPSE  o GY # 0  )  REWIND  LAPSE  . 

GO  TO  M'2«3«4'5»6*7«3«i0n«j<.  ..  ...  - . -  -  - 


si&n 


M  s  MTNQ:(;-N,LRL  <L>  J  . _ ....  -  ...  . - 


no  85  I  S  1,N  ■  . r'  ,:r’  \ 

Pack  space  l  . .  ’  .  . 


LBL  (i.  )■  =  LBUL5-N  ‘  '  “  .  '  .  '  ‘  «;/■< 

GO  TO  on .  .  ...  .  .  .  ; 


END  FILE  L  .  .  •  ' 

GO;  TO  00  .  . . 


LAPSE  *  -1  '  "  .  %  *• 

PRINT  49i  TAPUNlL) «  IDENT  f  U  .  .  . .  . 


REWIND  L  . .  .  .  .  «■ 

LALIH  = 


continue  . .  • 

TF  {  LAPSE.  EG;  0  )  GO  .TO  108 . -  .  .  .  ... 


PRINT  AS .  .  '  '•  *v* 

pAUSE  .  .  .  „  ...  .  .  ./  '  -  -- 


CALL  ATLAS  .  “  •  ~ 

GO  TO  \00  . 


CALL  STRIKE  . . .  >"  \  ' 

GO  TO  100  .  " 


CALL  DFFFND  .  '  - - ■' 

GO  TO  109  ' 


CALL  ASSESS 
GO  TO  109  .... 


CALL  SHELL  ...  - 

GO  TO  109  '  ,  .  .  .  ,  . . .  . 


CALL  MOTION  .  ..  -  .  . . 

GO  TO  109  . .  .... 


CALL  SURVEY.  .  ...  .  _•  - 

GO  TO  109.  ...  '  ...... 


CALL  COMAND  *  *  *  *  ;  ;  - 

WRITE  ( L B 9 409 )  NAME*  DATE i  JOB9  LBL  .  . 


GO  TO  101.  '  . . - 

FORMAT  (10 A 6 « 1 2 A 1 )  \  "  ..... 


FORMAT  t  lHl/////.////  . 

+  HX  39H  HEADQUARTERS  AIR  FORCE  SYSTFMS  COMMAND 

*  1SX  37H  ACADEMY  FOR  INTERSCIENCE  METHODOLOGY 

+  //  22X  23H  ANDREWS  A'lR  FORCE  BASE  •  .  .  :  ..  .  . 


«  ?6X  31H  MUSEUM  OF  SCIENCE  AND  INDUSTRY  . . 

*//  2?X  23H  WASHINGTON.  DiC.  20331  _ _  _ 


+  SOX  25H  CHICAGO' ILLINOIS  6063? 

*//////  S2X  JL5H  BEACON  PROGRAM  .  ..  .... 


+//  32X  5?H  BATTLE  EFFECTIVENESS  ASSESSMENT  OF  CONTINENTAL  DEFENSE. 

+//////  36X  SAG  //  54X  2A6  /1H1/) 


0 R u A T  ( 1 7r!  THIS  SECTION  IS  » A6»3Xt8A6i  10H  WITH  KEYS» 3  (4X^412)  tJY 
FORMAT  fSHOEND  OF  .  A692A^,.S0X,  i9H.  RECORD  POSITIONS  •  4 1 A  /  i  H 1 


FORMAT  ( 15H  mount  REEL  *0,A'6»UH  *  ON  UNIT  A6»13H  i  WITH:  RING  A 6 )  . 
FORMAT  (,?3H.  PRESS  START  .WHFM  READY.) 


FQRMAT(19H  SAVE  TAPE  ON  UNIT  A6«12H  WITH  LABEL  A6) 

END  _ _ _ ,  .  _ s . . 

c  IlkIUk  DTck  '  "  “  '  ““ :rr“r  • *•; 

BLOCK  DATA  .  . 

COMMON  /PAR/  LA  « L B  9  R A D 1 1 J S P  I E 9 RANDA i  RANDB e  ~  ~  '  ~  r_T^" 

♦  C.ONV9SCALAT9SCALON9POLAY.POLON9 _ .  .  ,  _ .  ...  , 

♦  R 0 IJ ND ,  EP S I ,  D I L  A T E 4 C P K  1 . C P K 2 » C PK3*  ~  .  ?  “  “  “  “ 

DATA  LAtLS9-RADIU£4PIE4RANnA<R.ANDB  /5t6t3437i7i46'8i3«:mS927v  . 

«■  4  0  5  74  6  5  0  t 9 . 7.4  359130/4  '  : .  "  ' 

t  CO  V9SCALAT9SCAL6NCPOLAT9POUON  /»4a-00^-lOQ>;-g^Or/<  .  .  ..  .  . 

♦  gjjjJND 4  EPSI , DI LATE « CPK UCPXi «TCPK3 ’  /2<>i 5 5341 f  9 « 5/ «  . .  > 


<\  *  „  *  A''  O 

•  V.  -  : 


■v  ^ 


0'  J?  *  J 


'  >■}.  *9%y*  i  r?,  ty/m 


SI BFYC  ATI.  TDK  DFCX 


SUBROUTINE  ATLAS'  '  “  "  ;  -  ..  '-'‘7  "  = 

COMMON  /MON/'  inENttS)  ,KFY(12)  ,LBL(4)  ,N(JM  .  -  -  > 


Common  vp  AD/  LA, I B  *  RADIUS,  DAT  (31?  CONV  iSCALAT  ,  SO’ALON  ?POLA'T « POLON. 
COMMON  /PAc/  NF,NAMT.YPl30i  ,MPZ(30*  ?VAU:(30)  ,RAnT(30)  ?  . _ 


+  W.'A’T  (T^OO)  ,WLON(lOOn)  ?VALU(1000)  ? SIZE  (10005  , SECT  U 000)  ,0l»AN(1000) 


i  ~ 

FQUTVAlpmCe  {NF-,FACI4>j 

nATA  nomp  mh  /?  degree  /5?*3  /  . 

< 

- 

* 

COMMENCE  '•  •  '  . .  . 7 . . 

IF{  MOONEY  (?),?)  .EQ.0  )  GO  TO  10  „ 

'  c- 

REA  .  ( V)  URL  (1 }? FACILE  ’ 

,  .  .  1.0  IF  (  NUM.LT.  100  )  NF  =  N.UM  .  . 

■  '*•  *"2  ^  ' 

• '  • 

'  r=  n  .  - . .  ■■  ~  . * 

FF  (  NF.NE.O  )  1=  MFZ  (NFL . 

'  *  \  ' 

'  ,  •  i  f  's 

•11  READ  ( L-A ? 41 )  NAME ? NV «NS?  INFMT,  I NOUT 
.  IF  (  NAME.EQ.MONF  ,)  _  GO  TO  ?0  . .  1 . 

.  y  .  \  ) 

"RITE  (LB?40)  .  . r  ' 

WRITE  MR? 4,1)  N AME ?MV. 5‘ N S.,7 NFMT^f  Mnu*  . 

<  v 

mf  s-fjhni. —  - .  ■  ' 

MAMryp  (NF)  a  NAME  . . .  . 

*  ■  * , 

Nn  -  1+1 

'  ”, '  ’  .  ‘  '»  t 

JN  =  NVRNS 

^  .  i  .. 

„  -  ,  ,,  ,  „  ,  „  t  .,« 

12  P;An  (LAi  INFMT) 

(VAR(J)  ?Jsl?JN).?KC  ~  ‘ 

c  1 

• 

T F  (  *C»MF.,0  )  GO 

TO  15 . - . 

I Ft  Y  EY( 6) « EQ  *0  ) 

WRIT?  (LB? INFMT)  t VAR  t J) ,J~1?  JN) 

*  .  v? 

„ 

J  =  0  .  . 

_  .  k  -  -<■-  „  ,  ...  , 

.  v  . _ 

'* 

- 

TO  13  JS  s  1,NS 

T  =  T+l. 

7  si  '  ^  0 

SECT  (I)  i 

13 

OHAM(T)  : 

TO  1? 

T  =  T  +  l 


VAL-UF  “  0, 
no  i  a  |  - 


WL  AT  {  t  )  z  SCAl.A  W  WL  A'.T  ( I)  -C0NV«AINT(  WL  AT  ( 1)1)  POL  AT 
•WLQ.M(T)  =  SC  A  LON<KW  L  0  N  U )  ~C  ON  V  #  A  I  N  T.,(  WL  ON  {  ?) )  )  *  POLON 


53 


>X'< 

«?*  ^ 

v  * 


V&LUF  ?  VALUE  +  VALUCir  . 


IF(  VALT{MF) ♦tQ.O.J  GO  TO  18 
writf  n-BtftOM 


'a  WLAT  { ME )  /RADIUS  " 

.  8  =  WL'OM.fNc  >  /RADIUS 


S 

S 


SIN (BIRTHED 
S IN ( A  V 


DO  }V  t  =  NR, ME  . 

A  a  WLAT (!) /RADIUS  ..  .. 


*  =  WLON { JO /RADIUS 

‘(.AT  =  cns(A)1  . 


SLAY  ’s  ST  N  { A  Y 
CION  =  COS  (R)  ... 


SI.  on  =  SIN(B)  . . . 

A  =.  CLONRCLAT 

• 

■  <? 

R  a  SLON^CLAj  '  " 

...  .  ST  7F  (  T  )  i  SHED'^SLAT  *»  THEDAS,  * 

* 

1  ’ 

SFCY(I)  =  THED»A  “  CH£n»SLAT  " 

•011AM  (T  ).  =-  CHfDoR  -■  SHEO^A 

• 

.. 

I F  (  VALY  (  NF)  .GT.-O.  )  GO  TO  19  ' 

THEn  •=  valu(t  ) /degree 

THED  =  SHEHtfSlAT 
WLAT  (  T )  .2  GHEOOSLON  . 


WLONTT)  =  "CnEDOGLON 
V/ALUM)  =  SHED«CLAT  . 


=  A 
=  B 


CLONRTHED 


‘  SLON<miED 


.  17 

•THEH  =  SLAT . . .  •  1  \  . 

WRITE  ( L  R 1 4  0  2 )  WLAT(I)  ?WLON(I)  «  VAfcU.t  M  VSIZE  i  !  >  « SECT  (IT  ,tf  QUANT!  I  - 

•  '  <* 

.18 

GO  TO  11 

VALT(MF)  =  VALUE 

' 

i  » 

TF(  KEY  ( 6’?- i ECU 2  )  GO  To  11 

I F  f  INOMT’tfl)  ,.EO»0  )  WRITE  (LB, 400) 

(SECT T  I  )  oQUAN  ( I  )  ,  Wb.ATTI )  ,  . 

l  ' 

*  *i 

c  WLON  ( T )  vVALU-UT  rSIZEM) ,  lsNB,NE) 

I F  (  vmohT(A)  iEba  j  WRITE  (LB«401) 

(WLATM)  ,WL0N(!)  .  'sNRoNE) 

20 

GO  TO  ] l 

IF (  WV( 2) « LT« 2  )  RETURN 

- 

IRLCD  ;=  LBLTD+I 

WRITE  (!  )  ‘URL  ( 11 , FACILE 

- 

40 

RETURN  . 

FORMAT  OUT) 

\ 

41 

42 

FORMAT  nXA6i'lX,'2I293A698l  i) 

FORMAT  (<5F8t  2) 

1 

43 

.  400 

FORMAT(/eXi2F12.2v8Xi-7ElO iW?~  "  ”  •  . 

FORMAT  (/52X  1-7 H  ATLAS  PARAMETERS  7/(1  XA3 AG c:3FT„.0«.F7.:iii  1TA 

»  V  7  '  ■  n  l  ^  •  >  ■  •  ,  r  -  V  7  —  p  •  (J  If  ?  ' 

-401  FORMAT  (/P6X  34H  LOGIC  COORDINATES  FOR  THESE  S I TES// ( 

40?.  FORMAT  (1?X,3FK  i  4  ,T?T,'3F  1  ? .  4 )  ~  ‘  '  - r 

403  FORMAT  (/SIX  1.8H  VECTOR  COMPONENTS  //V 

r“|N[)  -  rT-Tr-^-rr-rT-Tr-T - r-r~rrr - r~ - r - - - 


Mil: 


'  \s  * 

„*9  ; 


V  9\  V  >  vVt/yMA  X/  V  **,  '  o' 

'  ^  i  '=■'•" N  -i 

: > V' <  >  «V.  V*  <  '■ b-  ••  A&  ) 

4  &«•  st-  ..v •  * 

.-  ft  •  ckT  iakt 


-  .  -f  -  "  ;•  £  V  •'  A  ;  A",  ;<  r  1 


- 


SWTC  'ST0  7DS  DECK  •  '  . *'"*•*., 

~  stirtou-tTO"' stpjk e  ~  "  ;  '  **  ”  - v . *■? v” 

COMMON’  /MOM/.  ■tDENYr3Mxr,vT,RT4R.W^BYPASS-«KFY(r2j'UfJLU).  - 


Cfii^nN  /PA R /  L  A  i-LB.* R ADT-US «p I E i DAT  ( 7)  9 ROUNDS FP S I  ’  '  " 

Common  /rAC/  RF,^AMTYP>30).9MFZv(30)  ?VALt  (30)  «RAn.mO},,  - _  ' 


«  Y(inoo)  oXflO'-O)  «VAL  uonoj  iAPEA(lOOO)  iSECTUOOO)  «QUAMn000)  '  <■  77-' 

common /wwbvf  w4Mw7(i?)  .Ywj--.{2onj)  awHirfiono)  vPWmp.ooo)  iiPwdodoVl 


nrMrNSTOM  VAPK  aOOOj  ,PSP(lOOO)  VANPMOCO)  ,AMT  {aDOOI  rVAClOOO) 
n.iMrN‘S*ON  I-MpEY  ( 1?)  « NORn£R:(  IP)  *MWmi2}  ,VLA-(15)  «RStT(fi)  «PWHnKtl(&>  « 


+  A‘LY(,l?a^  '  ”  ^ 0 1«5 J  *WRA'DC12,15V «C0DECl3vl5)  «MWTGY{12«m  / 

_ ..  niM^stfv  ■  ,eEP(l2UVAL!M(15')rNRTd5)  •  WEAPON  (AO  13)  . 


'  TaUIYAiPN'uu  .  , .  ,M«IT) ■  «  (-VAUN'tRtQT *  .'JW< WEAPON}  . 

_ _ F3UtVAlFH.CE.  rxwHHoVSP'lQ  ,  (Y-WHt),P$P,ANM  9  (REQaMWTGTJ  «  tPWKoAND  .  . 


COMMENCE  '  '  . .  * 

TF  <  mopCkEv  (2) ,2>,£q,0  j  GO  TO  9G0 


R£An  oriBKJH'  VIFAPON  ' 

IF(  KEY (?)  iFQ#0  )  RETURN  .  .  ,  . * 


(•O'  TO  99% 

REAP  .  ALArAOO)  yN 


WR I  TP"  (LBrAOO)  YN 
RE (l.  A, 400,1.  NC 


WR-ITF  a^Aoorwc 
my  =  n 


3 

WRITE  ,M.R.,403).  NA 


y  unrx 


yndex- 

NORnpR 


NORDER 


9PABT 
NAME«KG9pA8T 


I,'l» 

no  0,11  kf  s  i,nf  ■  •  '  ■ 

If'  MAMF.EP.NAMtYP (KF)  }  GO  YO  912 

<  »  - 

"  911  CONY  INI  IF  '  •  '  7  " 

Go  TO  9t* 

>  1  r 

*  t  (  \,  | >  ■ 

912  MT  =  hit  +  1 

...  mrt(MT)  =  KF  .  .  . . 

-  %  '  ') 

.  >-  -  ...  *  -•  v  .v 

0 

DO  913  KW  =  1«1? 

RFQ^WVNTI  =  0«  .. 

WRAP :^W,MT) 
RE An  Yl A,A03) 

=  PABT(KW)  . 

N A  MF , KG  9 {RFLYYK W  9  NT ) 9  K Ws 1 , 1 p ) 

WRHE  (LBYA03) 
IFf  KG.-LE.0  ), 

NAME VKG 9  (RELY  (KW9 N T )  9 K-W?l „  1 2 ) 
GO  TO  910  . 

CO  TO  910 

9.14  ?1ZF  =  512»<T3  .  ,. 

' 

....... 

s  * 

READ’  (I.A9401)  CEP 

v  .... 

* 

..  W  -*  * 

1  '  “  *  '  " 

.  ,  WRJYP  (LB 9 401)  CEP  ...  ’ 

■ 

DO  953  KT  -  It  NT 

CODE ( 1 9KT)  =  1, 

■  . 

no  9?4  KW  s  1,1,2 

1  VGIJM  s  ABS(RFLY(KW,KT)  )  . 

- 

'' 

PO  9?5  J  =  hlO 
IF(ntG.CT.VSUM)  GO  TO  9?4 


'  BAG  =  GT.G  p  0  - -  - 

rnpr (Y-w-nuxt <  =  c onE-(KWoKT),  9  big  . 

YF<  C OOF ( 1 3  9  K  Y )  rOfTslf 7. E  )  GO  TO  923'" 
WRITE (LG*G?0)  KY  . 


.  v  .•>  ,  c  *>'  ; % 


%  ‘  %  4  d*  S*V.’'  ^  d 

.  >’  a  ,  yc.  *  i .  V"  *  ?  11  ;  v  r  f 

m  sv«r  7,  V'  77  '<■.  J3-  Vi  -'i 

’\\i  -•  '  • >v  -'A  ■-."77 

Wi,:.  .s-  *^vj*  5s . .  «w>  4^TPy.  >. 


923  CONTINUE 


mPNH  s  MPZ1NF) 


,<>v;  v 


Vat  IT)  » 

920 

A*-' I  (IT)  : 

, 

NW  =  0 

300  DO  399.  IV  9  U  1  ? 

- 

'  <t".. '  a 

’  «c.w  =  NO»DER  (LW) 

.  . 

"  r*  '  v  r 

A  '■  •  i 

• 

IF (  KwlLF.O  )  GO  to  940  . 

* 

n*  ’=  i^Aa0(N,w,KW) 

-  *  v 

>  -  1*.  «■> 

i  "** 

,  S 7  =  4, 9323-CEP.  *92  .  ’ 

^  '  ‘ 

MA,\ v  =  INDEX  (KW) 

IF  (  MANY ♦  Ef) .  0  )  GO  ?0  399 

.  -»>  •  ;  v*  6'  \ 

.76  %  1  =  1-.,  NT  ’  '  ““  “ 

‘ 

-  4/  ' 

KF  a  NRT(KT) 

* 

’  MPEG  «  1 

.  I F  (  KF;g.T.#i  )  MB  EG  “  Mp-Z(KF-l)  «  1 

*  1  ”  ,  <•  { 

mfNd  ?  MFZ(KF) 

* 

a  ' 

.  PRO  a  RFOs(KW#KT) 

(PRO)  329,313,312' 

AN  a  l+NEND-NPFG  .... 


s/*ul«F  =  APS(vALT(KF)  )  0 
VL6  (XT)  a  ALOG.  (VALUE’ 


<1»<*AINT(PRQ)~PRQ) 


aw  =  PIE  <>  WRAD(i(;W?KT)o92, 
PO  309  TV  =  MPEG?NEND 


v a p <  urr  -  n.o 

IF.(RFLV(vKW,KW»GE*Oi)  Go  TO  316 


REAP  (IA,40?1  KC. 


WRITE  (LB«*»02),’KC*(KSIT.{J)-,PWHDKU(J)-*Jal9-6) 
DO  964  J  =  1,6  . 


IF(<&-1T(JM  962, .964, 961. 

961  KE  =  MBFG-1+KS1T(J1  .  .  .  . 


VAPk  <KL'  a  PWHOK-U  (J>  '  • 

GO  TO  964  ...  *  .  ,  .  . 


962  KB  =  KE+1  '  ‘  ’ 

_ ,  =  NBEGr.l-KSITl  J) _ . _ .  .  ...  . 

no  Sb'3  K  =  KB,KE 

963  VAp’KCO.a  PWHOKL(J)  ...  ..  . 


964  CONTINUE  "  . 


W% 


RELYUW,KT)  =  “RELY(.KW,KTI  •  '  '  ' 

316.  PX  =  RELY  (KW.,KT.) -AI-Nt  (reLY  (KW,KT) ).  . 

V 

PX  a  AMAX1  (PX,,1.«PABT(KW)  )  ' 

-w,  -  1  ’  *  '  *1’  -  v  ,r,  ~  *  ' 

I r  (  RAnr ( KF)  )  317,3209315.  . . 

317  P<  =  RX  -  PXttEXP (-AW/SZ) 

GO  vn 

-  -  «  *  -  «•  *?  •*»* 

,  0 
■  <  A- 

315  AT  a  PIE'  *  RAPUKF).<N>2 . 

PK  =  PX  #  PKAV.tA.T ?AW,:52.)  _  •  _  .. 

..  _ _ 'L 

•  -  .  ~  r 

-  320  DO  ,-«26  IT  s  MBE69NENO 

T  ~  1  i  V  .* 

.  !FtVA(TT)..CT. VALUE)  Go  TD  322 

. 

’  VAPK (IT)  £  0,  '  ------ . 

GO  YD  328 

| 

322  IF(RADT  (KF)  .NE.O.)  GO  ^0  .326 

321  AT  =  A R E'A  ( I Y ) 

' .  \  *'*  ■-  -  .  ‘1 

PSP.f- r.T)  ,  S  1  ,.0»P-K«  ( I  .  “VAPK  «m) 


VAPK  (IT)  3  VA(n)9{l.*?SSP(TT) ) 

2-7 


°5P(tT)  s  *-5f*(IT)'*»fl..yAREA;'tT)  > 

VAPf  t ! T)  =  VAPK(rT)/4RFA(?v) _ 

326  C.ON-T  I N*if 

329  CONT ?Nmf 

330  VDMAT=  0.0 

no  339  KT  S  1  ,MT. 

«F  a  MRT(kT) 

r  1 _ _ 

"  TH  rTn?. :  >  Nr)£6  =  ^fz-ckf-u  +  i 

_  ,Mrvn  -  Mr?(KF> _ _____ _ 

^  r  ( 9  ww  i  K  W » K.T )  )  339 , ...  JO ,  334 
334'  nO  V3A  TT  s  MREGiNEND  . 

TF'tVPMvi -VAPK'UT)1)  336,338,336 
336  VPM4X  =  V.APK  (IT)  , 

•  jT  v  ~  “~" 

LT  =  KT  _ 

333  COMTINWF  :  r^~r 

339  COMTIhtUF 

!FUVBM4*)  350*350,340  '  '  '  " 

3r40  VAPv'(jT)  a>*0 

A*N"  ?  ALOGtvA  (JT)  )**VLG(LY) 

342  P$  a  PSP  { JT)  ‘  ........ . 

'  RN  “s  AKNP/AMA'X'l  ( *  001  «**AL'0G  (PS)  ) 

IP  (j^Fy  (5)  ,f,T »  Q:)  ROUND  a  RANDOM (ROUND). 

M  =  INTfRN  +  ROUND}' 

..  SWIFT  =  CdnE(KWiLT)  _ !_ 

-“L'PAS't  '=  1  IT  (RciQ  (KWVLT)O  “  “ 

_ N  =  -m-a.XO('M,  LEAST:) _ , 

MOST  =  TNT.(RELY(K-W#LT)  ) 

_ IF(  M.LF.Hnsr  )  „GQ  TO  34 1 . . 

'  KF  =  NR V  {£ T) 

WRITE  gRj.335)  NAMTYP  ■(*£)  ,QUflN-tJT)  ,N,KW 
‘  "  H"=  MOST  ' 

347  ,  N  s.  MI  NO  (NjMANY),  . _ „ _ . 

"  ‘  VA  f  JT)  s  VA(jmpS<K>N 

IT  JT;)  a  AMI  (  JT)  *<•  FLOAT  ~(NH>  SK!FT 
many  =  MAKJY«N 
IF  (MANY)  39.9,3999  330, 

150  /PRO  s  MANY  ‘  '  '  :  “ 

„  .  7* NT  liA _ _ _  .  . 

'  '  MP  a  0 

351  VP.M-AX  a  0,0  _ __ _ ..  „ 

"  DO’  359  K 7  a*  1 , N  T  . 

„ KF  a  MPy/KT) . 

NRP.5  ~  1  '  . 

TFC  KF,GT,1  )  NBEG  =.  MEZ.HCFwir  »  1 
»FMn  s  MF7,OCF> 

IF  (RFQ  (VVhKT)  )  359,354,353 _ _ 

35  3  R F Q;f K W ,  K  T )  s  -1,  '  ‘ 

r  -'.  w  a  ^ 

,  1YU  I  *'•  ...  . 

154  no  356  ) r  a  wnFG9NFNn  “  '  ' 

IF  (VPMAXrVAPK  HT)  )  356,356,355  . 

356  V PM 4X'  -  VAPK  (lfV)  " 

_ JT  a  TT _ _ _ . 

356  CON  T?  Ml  IF  . 

359  roMVlNUF _ _ _ ■ 

!fT.  PMAXJ  3 TO  <.  3 7.0 ,300  . 

AKNP  a  »A10G(VPMAX-) _ , 

7FC-AKhP'»  ZANY.  -ZPRO)  361«3TL*371 
.  ,  .  .  2-8 


3.60 


M?  *  \ 


VAPV(JT)  “  «AKN?~2fi.O 

A  M  -  1.0/  /.  maxi  ( <  r:, ~;.i  os  t  o  -to 


f  (Mpj  5-/0,390,371 

'  '/  ,\:  £  —  7  ;>  n  n  /  7  A  M  f  - 


*  5  )  N 

MP7  CCF 


i  *  (  1  i.  Q  (  '  n  <  K  T  )  j  t  o  0  o  )  00  •  (  Ik  V 

On  37  9  TV  «  H3F0  -  MFMT1 


I c  {  VA*K  C I T 5  «0f  »0C  5  GO  TO  379 

IF  !  VPmAV.GF.VAoXMT)  5  00  Yn  979 


IT 


V*».M  ax  =  V  A  P  K  (7  T) 

?  7  9 

OMTTN.17 

„ 

369 

i 

f 

CONT  T  \'t!f 

IF C  VP v AX  «  FO  0  ~40 . 0  )  GO  YO  300 

V  =  (AKNP+VAPK<jT)+20.)  <>  ANP(JT) 
•/Are  <  JT a  =  n. o 


TMT(EN*t.O> 

r:  T  N  Y  ( R  p  I  V  f  K  .  L  Y )  \ 


IFC  M.LF.MOST  >  GO  TO  .366 

" F  r  NRT H  Y1 


WRITF  (LR«a351  NAHYYP(KF)oQUAN(JT)9N,KW 


Ml  NO  (N 9 MANY) 


* 

SHIFT  =  cnnE(KW,|.TJ 

V  A  { J  T )  =  VAC  rf)  0  F”P  {-AN/ANP  C  JT)  ) 

1 

4 

AMJCJJ)  =  ANKJT)  *  ANtfSHIF? 

MANY  =  MANY~N 

- 

i, 

1 

IF  (MANY;  399,399*372 

330  V R  I  7 F  ( i  *3,333)  M A M Y « K W 

j  GO  "n  3S0 

1  590  WRIYF  C  l  B  «  A  3  9  $  M  A  M  Y  «  X W 

V  4  L 1 '  fc  =  0, 

V  5UM  — 


no  941  XT  =  ] ,  NY 

XF  -  MR ” { e  7 ) 

• 

— 1 

I  MBER  =  1 

I  IF  C  i' F » G ? » l  )  NPEG  s  Hf/  (KF“1)  +  1. 

- 

m.  Nn  =  M’r  /  CKF) 

j  WRITE  n.P,342)  mAMTYPCKF) 

“ 

.  ,  j.. 

j  WRITE  (LH,890)  C  Vu  U  Y)  *  j  T=Nb£G,NY:NO) 

' 

. 

" 

' 

1  r a 3  vi  g  urn  =  o.o 

I 

V  A  L I M ( K  T )  =  0, 

~ 

00  GA4  IT  s  NRFGoNFMD 

1 

VIGH'T)  =  VIGCKY)  ->  VAf?Y) 


Z  - 


’*'!**■?  til 


/SHIP 


.  ■ '  j  -  •«  mi  n  -flo* y  imuuvi)  j  *skikt 

VO  Ms  M+M, 


MW "i  u  T  { s  v  is *f )  vMW |  Y  ( K W } 


<:'lv  -  L V »- 1 

r p '  v'  -  t  ->  ■c-  :v»nFv 


*  -  f"G  -  (V  ;  J  H  /Ra3IU.S) 
';T'*  - 


■  XI  IT)  -  C J£V/EX 

r  '  C;  r<  -  YCTJ  -  CDF V 


v  fa ;  j  ~  uly 

?a2  i t  i  n.i y,r,:E  )  Gn  t"  gay 


■  r  !  *.<T"UWJ  ,LL  ,0  )  GO  TO 

Ir  rEY(6)  .®T.OJ  GO  TG  5&6  .  .  .. - 


1 1> a  -  iNnpA»^w)  ->  i 

’^r**r  (vWHPino  ,XWHD(IWJ  flW=!WF-i!WE) 


/  h  u 


M  */  I  1  *  * 

/,  n  ,<  ,/v  , 

•iri  U  IHDF  f •'?’/) 


:.i£  s  IWB  v  mWIY.'KH) 
GO  Y'"  G  A? 


__"_  XWHnnWnj 
YWHnuWPf 

_ IF(STTFjGE,A 

STTF  =  SITE  +1.0 
548  YT  =  YT+PEV 

X T  =  xT+XPEV 

GO  TO  849  .  .  _ „ 

'537  X  T  s  XTl  i  ) 

YT  =  Y  (  I  T )  _ ' _ , 

T  F  fp  A  D  T  (XT )  .  roTrf.  )RT  =  SORT  ("AREA  (IT) /PIE")  t 

534  rso  536  KW  _=  1  ,_NW _ _ ; _ . _ 

JFf>  WIT(KW)  f~536 » 536+  535 

_5  3  5  _RW=  WRAP (KWiKT) _ „ 

Cali.  coord <mwit(kw)  .index  <Kwn 

536  CONTINIIF 

589  Continue 

_ WRI_TF  ( L P^8 58 )  (MWTGT  (KW.KT)  »KW=1  »NW) _ . _ 

00  590  KW  s  1,'NW 

590  MWZ(KW)  =  MW7 (KW)  +  MWTGT (KW.KT) 

Wconttniif 

WRT-TF  (LBi897)  :  * 

HO  997  XT  s  1 i NT 

_ .  X  F  -  MRT  C  K  TJ _ _ _ . _ ; _ _ 

997  WRIT?  (l. P>853)  NAMTYP  (KF)  .KF?  (MWTGT  (KW.KT)  .KW~1 ,NW) 

WRITE  (LBi,859)  MWZ _ . _ , 

MEMfl  =  0 

po  XW  =  1 ,NW 

MFNP  s  NpND  v  MWZ  (KW)- 
.993  MWZ  (KW)  =  NENP 

998  REAP  ( L A 9 40-1 )  PABT 

WRITE  (L Rj»4 0 1)  PART _ .  _  ; _ _ 

REAP  (LA t 401)  CEP 
WRITE.  (LB, 401 )  CEP 
I.  =  0 

NENP  =0  . 

Do  995  XW  =  1 i MW  ’ 

MB E 0-  =  MEMD  +  1  _ .  , 

NFMH  =  MWZ (KW) 

IF  (MEMP.LT  iNDFG)  GO  TO  998  . _ ; _ 

I F  (  KEY  (7) : EO, 0  )  T  *  INDEX (KW) -NEND 

_ S?  ~  0.84933OCEP  (KW) _ „ _ .  _  . 

ZABT  =  0. 

DO  994  j  =  NHEG.NEND _ ,  _ . _ 

CALL  RANORM(DX.DY) 

JI  =  J  +  T _ - _ ; _ _  . 

EX  =  COS(YWHOUI) /RADIUS) 

_ XWHD(J)  =  XWHP(JT)  +‘S7^DX/EX _ .  . _ 

YWHn(J)  =  YWHD(JT)  +  S7.«DY 

_ PWK(J).  =  PART  (KW)  , _ . . . 

IF(KEY(5>  .GE.2)  PWK(J)  =  A I  NT  (PWK  { J) +RANDOM  (ROUND)  )■ 

994  ZABT  =  ZABT  +  PWK(J)  _ .  _ _ 

TmTP'Y  ( 6 )  ♦  G  i  «  i" )  GO  TO  99D  ' 

WR I TE  (18,890)  K W. ZABT 9 ( YWHP( J) ,XWHD ( J) 9 PWK (J) . JsNfigO . NENP) 

995  COMUUMf  * 

_ I F  (  K-EY  (?)  ;LT.  2  )  RETURN _ . _ _ _ „ _ , 

LRUS)  =  LRL  ( 3)  +1 
WRITE  (3)  LRU 3)  «  WEAPON- 

PETURN  •  *  *  ~  ~  ” 

400  FORMAT (8H  CARDL8L < 1216) _  ; 

AOTTORMAnfiH  CARDLBLol2F6»2)  — ■  — 

2-11 ,  .  '  . 


’  o  < 

V  -  C 


4  Y  -  * 

A  "s  4-.*S  'V»  '  W 

'V  0  O  V-  V' 

\\  :  ■Vyfjt'?5 


403  FORMATOH  »A6tA‘l*12F6.2) 

020  roRMATt/ /32H  CODING  OVErFI  OW  ON  f ARGFT  T'YPF  I  31 _ ^ _ 

FITTo"^TT ?XAFT6H  TARGFT.i4XA6iI0H  RE0UESTEDI4.8H  WHDTYPFU) 

A3  A  FORMAT  (11HQTHERF  WFREi6,3X16HWARHEAPS  OF  TYPE  I  3«  3X  3HLFFT  .,'»//) 
83rT0RMAT(irH0THfRF  AREl6,3xl6HWARHEADS  OF  TYPEI3*3X0HLFFT  .  /•/) 

030  c0RMAT<  1  HI  30 X  4QH  f  HE  EXPECTED  VALUE  OF  TARGETS  LFFT  APF  /) 

m2  format  (/^xag*  oh  targets  /)  . . . 

*40  FORMATS  ,•  X«  10(F10,?)  )  >  . 

84l  Format  {//SH  50m  FTHTTgR  AFfEft  t37X10H  PER  TYPE  /(1AX»10(F10.2)  )  ) 
043  FORMAT ( //5H  SUM  Fli.?»8H  5F(f>0RE»37X10H  PER  TYPE  / U AX , 1 0 (Fl 0 . ? ) ) > 
A^O'TORMAT (lHl  37X  44HTHE  ~Rf3ULTS  CF  THE  TARGETING  ARF  AS  FOLLOWS  /) 
052  FORMAT f//2XA6«8H  TARGFTs(42X25HNOK*ER  OF  WFAPONS  OF  TYPF  /  POX  ♦ 

+  rHT7xTHt7X  IH37X lH47XlH57X 1H67X IHTTXI H87X 1H96X2H1 06X2H116X2H12//) 
053  FORMAT(PXAA,2XI4i7X12I0)  ... 


050  F0RMATI//22H  TOTAL  WEAPONS  TO  TYPE  I77TM8/) 

097  FORMAT  (1H135X50H  NUNRFR  OF  WARHEADS  BY  TARGET.  TYP'F  AND  WFAPON  TYPE 
♦  ///?  1  H  TARGET  7  WEAPON  TYPE  ‘  " 

>7XlH17XiH27XlH37XiH47XlH57XiH67XlH77XlH87XlH96X?H106XgHn6X2H12//> 
859  FORMAT  (//^4H  THE  GRAND  'TOTAL  WEAPONS  I5*U!ft/lHl//) 

861  FORMAT  (6X«2F9«2>  > 


SIBFTC  MOVKDK  DECK  - - - - 

_ SUBROUTINE  MQflnN 

COMMON  /par/  la, lb, RADIUS, PIE  ' — - - — " 

- COMMON  /fa c/  NF«NAMTYP(-30^  ,MFZ(30>  »VALT(30t  ,RAnT(30i  . 

+  Y  (1000)  » X  ( 1000)  ♦VAL(lOnO)  ,4(1000)  id  (1000)  .CUOOO)  7 

- COMMON  /WHn/  MWfMW7(l?) ;Vw(2000> iXW (2000) *PWK(20fld) . IPW'RnOM 

COMMON  /MAli/  LPZC12) - - - - 

- _ DIMENSION  NRP(l2)  ,DSLtM(l?)  jMANV (30)  it>WELT(30)  ,L(100n) 

WlVALFNCE  (LiB)  - — - - - 

COMMENCE 

F  d'EAn  (OT5Ti 

_ DO  2  KM  s  1,NF 

f"F  (  mamTYP  (KM)  #FQ#NAME3H  )  00  TO  3  '  - - - 

2  CONTINUE 

“  RE  TURN  ; - • - 

3  Write  (LB«41)  NAMESH^KCiNRP 

TF(  KC«LT#3  j  GO  TO  9  - - ~ ~ 

CLASSIFY  FACILITIES  By  SECTOR 


DO  «5  JW  =  1,12 
DO  A  KFs  J,NF 
TP  r  NAMTYP  UP)  TfBTRrF  ? 
continue 
"go  to  5 

IFP  s  1  - _ 

IF(KFiGT.l)  TFB  s  MFZd 

I~E  s  -MF-?  (KFl _ 

T  =  1,30 
MANY ( I )  .  =  0 
DO  3D  I  s  IFR,IFE 

_  Lin  as  locatf  ( y ( n  T 

II.  =  L(I) 

MANY  (IL)  =  MANY(TL) 
HTrTTF  (!.BiAo4)  nAMtyT? 
WRITE  (LBi^Ofi)  (C  < I J ,L 
“Continue  - 

GO  TO  ) _ 

TFT  K  C  •  F 0  •  2  )  GO  TO  8 
REAP,  (LA, 42)  DSLIM 
WRITE  (Lfl,42$  DSLIM 

IWE  a  0 _ 

"DO  10  KW  a  j  7NW 
TWB  =  I  WE  «■  1 


MFZ(KF»1)+1 


71  CONTINUE 
GO  TO  tO 


2 

I F  < 


t  PE  =  MFZ(KP) 
no  16  TP.=  1*30 


nwKiT(.lP)  =  n. 

18  _  MANY(JP)  SO 


Ft  XC.FO.O  )  GO  TO  1 
COMPUTE, LAUNCH  POSITIONS 


Ft  n  s  I,  TMJKW)-*  EG « O*  )  DSUM(KW)  a  10800. 
nSMTM  =  10800. 


^0  12  TP  =  IPfHIPE  . 

ns'  =  nisT(Yw(iw)*xw(iwi,Ynp><x(iPM 


I 

Ft  ns, GE.DSMIN 
HSMTN  s  ns 

)  GO  To  12 

.  . 

12  C 

T p W { I W )  s  IP  < 
•ONTT  nme 

-  IP8  ♦  1 

IF(  nsHIN.GT.nSLlM(XW)  )  PWK  f  IW)  =  PWRUW)  ♦  10i 
!  GO  TO  15 

IP  =  lIPWMW)  +1 

11  MANY  ( VP )  S  MANY  t  IP)  4>‘  1 . ' 

LPZ(KW)  =  IPB  «•  1 

" 

WRITE  (LB.40U  NAMTYP(Kp)  «KW*MANY 

CROSSINGS  OF  BARRIER 

16  DO  21  IARC  s  IPp ♦ IPE 


PWK  UW)  ,GT.O. 

JW  =  I P W ( T W )  «■  LP2(KWV 

xi  =  cos  (xi  jw>  fr « JwJ  /Radius) 

VI  =  SIN  (XUW)  /RADIUS)  OC0S  C Y ( Jw>  /RADIUS) 


NtY(JW) /RADIUS 
1  =  COS  t  XW ( IW)  /RADIUS) #C05 (YW  UW)  /RADIUS) 


I'WmiCKCKUYiifJ  11171:71  ill 


Z2  =  SIN  (YW  (IW) /RADIUS) 


m  =  71«X?-Z2«XI  -  - 


CC  =  X1PY2-XPPY1 
.  ,  .  JARC=I«F  „ 

•  '  \ 

AHI.Lf.SsJlLSXLY7t'JW.h«  XfJWL.-YWXI 
,  ,  .  n020TAoC=IPR, IPF 

WJ-..XW.I-I.W 

)l  ----- 

T>(Y(  JARC)  »X’{»fAKC>  VY MARC)  tX  (I ARC  1TI 
BP<>C  1 I  ARC) »CC^3 fTARC)  .  . 


CCoA(IARC)-AA#CfIARC) 
AA#8 ( I  ARC)»BBoA  f IaRC) 


J**  -T-t'  -Cv"*}?1**.  ^  ' 


r>0  304  j  .s  1,2 . .  „ _ 

_  f  EM»sp'f  S  T  J  A  L  A  T  ,  A  l  ON  G  *  Y  Q  A  R  C >  »  X  M  AR  C )  ) 

•TCMPianisr  t ALAT, ALONG, VtJARC)  ,X(;JA3C)) 
t  F  ( TEMP+TEMP 1 ,  GT,  1  •♦BMll.ES)  GOT0303 

.TEMPsHTST  ( ALAT, AL0NG,V  <  J W ’)  ,X(JW)) _ _ _ 

~  tEMP  1=I1I5T  (ALAT, ALONG,  Vw  MW)  ,XW(IW) ) 

TF  (TFMP+TEMPI  •  GTM.«+fiMlLES)G0TQ303 _ __ 

CROSS=TFMPl/nSLTM{KW) 

XLATsAT^T  (ALAT/60* ) ♦(•Ol  )PAHOD(ALAT,60») _ _ 

r^L0N«=AINT(ALONG/60.)+(.01)*AMO0<ALONG,60.) 

WW=TW-IWB+1 _ • 

IAR=TAPC-IpB+l 

URlTEtLP,402)kW,IWW,IAR«CR0SS,XLAT<Xt.ON6 _ - 

MamY(IAP)=MANY(TAR)-»-1 

_  nWFLY(TAR)  s  DWELT  MAR)  ♦  CROSS*  1 1  ,-PWK  MW) ) 

V  A  LMARC)  s  VAL  M  ARC)  ♦  1*  -  PWK(IW) 

GO  TO  25 _ _ _ . _ _ _ 

303  ALAT=-ALAT  :  ‘  ' 

AL0MG=-St0N  M 0300 «"ABS (ALONG) , ALONG) _ 

304  CONTINUE 

JARC=1ARC  .  .  .  .  , 

20  'Continue 

25' continue  ...  . 

TP9E"IPE-IPB*1 
no  ??  T  AP  =  1  *  I PBE. 

MARC  =  TAR  IPP  -  l 

22  HWElTUAR)  s  nWELTM'Afn/VAUIARCJ 

— wrrrrLTr^TSNmYis77pi-i((w, dwelt - ~~ 

WRITE  (LB»409)  IP0E,MANy 


WRITE  ( L  B  »  4  07 )  (VALMARC)  •  I  ARC*IPB»  IPE) 
10  CONTINUE 


GO  TO  1 

41  F0P«AT(1H  »A6*H  ,12A6> 


42  FORMAT (AX, 12F6.0) 

404  .FORMAT  M1H0NUMBER  OF  ag;9h.in  MESH  A6  / (3 ( 1 0X516) ) ) 


401  FORM 


Hif.iffLii7ra^y3:nu>ni  T.TTiTiiirff.u.rin 


4,75X16/(3(10X516) 


402  FORMAT ( ah  WHDTYpEI3,5H  PEN  I3*12H  CROSSED  ARCI3*10H  AT  MINUS  F6, 


>10H  mINS,  AT  2FA,2) 

403  FORMAT (?5H00WELL  TIMES  FOR  BARRIER  A6*17H  WEAPONS  Of  TYPE  14/ 


♦(2(inx,5F10,l) ) )  t 

407  FORMAT (/(3(10X5F6»1> ) ) 


403  F0RMATf/U0(2XA6,I4))) 

409  FORMAT M 14X16/(3(10X516) ) )  .  .  * 


♦  ’ 


SU8RCUT  INE  DEFEND 

*  CCWCN  /MON/  I  DENT  ( 8 )  »K£Y  ( 12 ) , LBL I 4  ) 

CCWCN  /PAR/  LA»  LB» RADIUS*  _ 

CCFNCN  /FAC/  NF,NAMTYP(30)»MFZ(3'2)»VALT(30)  »RADT430)»  "  . . 

♦  Y (1000 ) »X ( 1000  )  *  F(  1000  ), AREA ( 1000 )  , SECT (1000)  ,  GUAM  ( 1000.) 

-'CCFf^CN  /WHO/  NW,MWZ(  12 )  ,Yw(2000  )  tXHiZ0WY7PtiKfZ60B)~tlP\F{20&Z  J 
DIf.cNS  ION  R ENTRY!  12) , AZ ( 12)  ?  FBOGY l  12 )» DWELT (12)  »PDETI12)  ,PKCHU2> » 

- +  COVER!  12 ) » DR  I FTX ( 12) »DRlFTY( 12 )  ,  WTfcD  UTfTZ  WTNCTiT)TWK'n2  )YZ  UK  Q  2 }' 

+  » WEAPON ( 8013 ) » LSEC ( 1000 )  ' 

""  ECU  I  VALENCE  (NW,  HEAP  INTVTrSFC,  SECT!  '  *  ' 

COMMENCE „  -  ■ 


REAC 
WRITE 
'REAC 
WRITE 
READ/ 
WRITE 
'READ  ‘ 
WRITE 
DC  5 


(LA, 43) 
(LB, 43) 
(LA, 42) 
(LB, 42) 
(LA, 42) 
(LB, 42) 
(LA, 42) 
(LB, 42) 
KW  «  1.1 


RA  IDT 

RAIDT 

RENTRY' 

RENTRY 

AZ 

AZ 

FBOGY  " 
FBOGY 


DC  5  KW  =  1  vNW  '  .  2.0 

CRIFT  *  l./TAN(RENTRY(«W)/57.3)  21 

CRIFTX(KW)  *-DR I FT*S IN ( AZ ( KW ) /57 .3)  :  '  "  22' 

5  CRJFTY(KW)  *-DRIFT*COS(AZ(KW)/57.3)  23 

TCTREAC  ‘  ( LA, 41  )NAME,$70im'©* S AT U RVCKA NSTA'VATUT TR K r,T( LIH7 DISCS PVT-iTGW  24“ 


DC  11  KD  =  1  ,NF 
I F (  NAM6.EQ.NAMTYPOC07 
GC  TO  90 

WRITE  (LB,41)NAMEr$',0'E 
REAC  (LA, 42)  COVER 
'WRITE'  (18,42)'  COVER 
READ  (LA, 43)  PKQH 
WRITE  (LB, 43)  PKQH  r  " 

REAC  (LA, 43)  PDET 

WRITE  (LB, 43)  PDET  “  “ 

REAC  (LA, 43)  DWELT 

WRITE  (LB, 43)  DWELT  ~~ 

DC  15  KW  =  1,12 
‘  COVER (KW)  *  COVERTKW)#*2 
WK(KW)  ■  0. 

'  WIND(KW)  «  0."* 

ZWINO(KW)  *  0. 

“ZWKIKW)  *  0i~  ' 

IDB  *  1 

IF(KC.GT.l)  IDB  w  HFZCXD-TT"*-!' 
ICE  *  HFZ(KD) 

DC  22  ID  *  IDB, IDE  '  r - 

IF(F(TO).LE.0.5)  GO  TO  20  ■_  . 
KVJS  =  0.  1  ~ 

BOGIES  *  0. 

"  IWE  *  0 

EX  *  CCS ( Y (  ID l/RADIUS 0 
'  AVDWLT  *  0.  “ 

P 0  30  KW  «  1»NW 

- W I  NO  C  K  W  )""*~0T - - - — — 

I WB  «  IKE4-1 


_GO~TO 


■OH  C 


l.WE  =  MWZIKW)  • 

If  UWE.LT.  IWB)  GO  TO  30 
7F-(-CCV.ER  ( KW  l  .LE  .0  • I  GO  TO  30 

XO  =  X(  ID)  +  DR  I FTX  ( KW  )  /,EX».H  I.GH  _ - 

YC  =  VHD.)  +  DRIFTY(KW,)*HIGH 

DC  50  IW  =  IWBtlWE  _ . 

If  (PWMIW )  .GE.  1.  ).  GO  TO  50 

RANGE  =  (  (XWMW)-XD)*EX)*»,2  4  (YWUWHVP)«»2 
I.F(, RANGE. GT.COVER(KW) )  GO  TO.  50 

IF(  S. 1,6.0.)  GO  TO  51  _ _ 

Tf(  IPWUW).NE.LSEC(ID)  )  GO  TO  50  '“““1  : 

WlNO(KW)  *  WiNOiKW )  ♦  1.  -  PWKUW) 

PWKUW)  =  PWKUW)  -  1. 

CONTINUE 

I F {  W IND(KW) iEO.0. )  GO  TO  30  ~  "  ~  " 

ZWIND(KW)  *  ZWIND(KW)  4  WIND(KW) 

PCPRES  *  0W6LT(KH)/RAI0T 
■  W PRES  *  WIND(KW)«PCPRES 

WOBS  *  WPR6S *POET (KW )  5  T 

WVIS  =  ,WVIS  +  WOBS 

AVOWLT  *  AVOWLT  4  WOBS*DWELT ( KW)  “  '  “ 

BCGIES  *  BOGIES  4  WOBS  4  FBOGYIKW)*U.-D1SCRM)*WOB5 
CONTINUE 

HRQT  a  0.  ’  • 

IF'IBCGIES.EQ.D.  )  GO  TO  25  ~ 

PTQD  =  EXP(“0.5»CBOGTES/SATUR)**2) 

&T  *  eOGIES*PTQO 

AVOWLT  =  AVOWLT  /.WV I S  _  _  ■ 

HRQT  *  HUM#TRKT/AVOWLT  " 

ZGCT  =  AMIN1 ( CHANS*  BT*HRQT)*RAIDT/TRKf 

ZFH  =  ABS(AVAIU*F(IO)  r  ~ 

ZH  *  AMINl(ZGCI.ZFH) 

I F  l  AV  A I L  .  LT .  0 . )  F  (I D )  »  AM AXl t 0*7HTD)-ZH )  ‘ 

HRQT  *  ZH*TRKT/(RAIOT*Bt) 

I  WE  ~  0  *  . 

DC  60  ,KW  =  1 1 MW  . 

' WK(KW)  =  0.  ^  ~  ~  ™  '  "  '  ”  ^ 

IWB  •*  IWE  *41 

IwE  *  MWZ(KW)  .  “  '  "•  " 

.IF (  WINOIKW). EQ. 0.)  GO  TO  60 

XC  *  xnbj  4  DRIFTX(KW)/EX*HIGM  ;  ““ 

YC  *  YMO)  +  DRIPTY1KW)*HIGH 

PT  *  PTQD«PDET(KW)  *  "  V" 

HCT  -  HRQT«0WELT(KW)7TRKT* 

DC  70  IW  *  IWBi  IWE  *’ 

IfiPWM  IWI.GE.0.)  GO  TO  70 

RANGE  *  HXW(  IW )-XD)«EX )**2  4  ( YWUWI-Y0I**2  ■ 

PS  *  1.  -  PKQHiKW)*U.~OEGRAD*RANOE/COVERIKV)) 
wkb  *  IV 4  pwkuwT"  .  : 

PWKUW)  *  1-  *  ?H.*1  I W 1  *~r 3»»Hv7 i  s 

I F I K EY { 5 ) . EQU  >  PWK  ( I W )  -  Al  NTTP  WKil  W JHANCDMUiTT - 

WK(KW)  *  WK(KW)  4  PWKUW)  -  WKB  ... 

CONTINUE  ~  “  “  ;  7 

ZWklKW)  *  ZWK(KW)  4  WKiKW)  ' 

~CCNT INUE  ~  '  ‘  ‘  1  ~  T“ 

JD  *  ID  -  IDB  4  l _ * _ _ _  ■  ■  ,  '  ■  „ 

2-17  .'j  *.  *  ‘  v;  *? 


.56  1 
57  1 

50.1 

59 

60  | 

6.1  i 
'62  l 

63  7 
,64  I 

■g:i 


70  if 


.87  1 

'  88“  " 

89 

'  90~;  ■' 

91  7 
"“92  “1 
93 

'  ‘9^  «, 

95  ii: 
96'  iS 
97  v 
'  98-15 
.99  g  < 

1001' : 
101 

102.  % 
10.3  ;4 
".104“' L 
105  .C 
i0,6~, 
107 

'108“t  i 

m  m 
'■not 

in  if 


*,  •.  */  -I  v'r'  A 

'  1’  ‘  76  ‘■-7  <  '0-  4 

'  '  •  O'  "o 

.*  c  ^1““-  '*  » 1 
;4  '  -p  -  '•  Jri  '  O'- <>'  l,  „  .  *-  »  J  -  ' 

>1 4,;17.  A-.j.  f  -If  A  -.I 


S  IBFTC*"A'S‘WR'~QF'C'k"''  '  '  - - - - ” 

. SUBROUTINE  ASSESS  .  ...  .  ..  .  . _ ^ 

common  /mon/  TWtT«TTkFy  2  i' Tfc&T (4 )  ■ 

frlw  LA,,UB,RAnlUS,RtE  ..  .  ...  _ . 


NF»NAMTVP  (30)  »MFZ(30)  »VAUT;130J  ,RAnT(30) « 

♦  y  (icon)  ,x  aoor  )«vAnooo)  «areauooo)  ,Sf<cThooo)  ♦oiian  m  don) 


0) ,PWK(2000) «!PW(2000) 
DIHF.MSION  PS  I  Mfl) ,CAL  !  ift)  «,REN  1 14)  «PSIR{  i8«3)  ♦CALR  1 1 8  *  3 )  »UPW  (13  /  ♦ 


PSIOU8 

3),PK<3 


ST  C  AOOliOOSCAK  400i:,DOSREN  (  400)  »FACILE (6121 ) 
DIMENSION  STRINGtlO)  ,ian(A«3)  ,MAP(l200),PF(9i  iPV(9)  vTTTLEO) 


12) «Y!E^n(12) ♦niRTY(i2f  «HEAT(l2) , WORD (6) «LETTER(fi) 

EQUIVALENCE  (NF< FACILE) « (WORD, LETTER)  _ , 

EQUIVALENCE  (DOSPSI, MAP} ,(D05CALiMAP (401) ) , (DOSREN,MAP (801) > 

.  -DAT A  TITL  E  ( 1 )  /18HSURVI VrLAST  FALOUT/ 


*  -  >  i  ni 

3«  t4  •  *  S« ,9*  ,13 
■i..  CAL/O,  i2»  ,  4» 


+  7333* ,8667. , lOOOOi/ 

_J)ATA  PSTR/ 1 4 « 4 9 i  r7.i  %6i  lt!>i3«4»0*3*5«2«9*J?«S*2f'3<i2i  ,1,6,1*2«;7,,6« 


+  .5*..24  0.»22.K14.8«.tOt,7»Q.2»T,»5;t4.8»,3v9.3;i»2;fli2;3»i;i«AR0;* 
•*•■  42;  ,36.,26o8,l6«6,l2:*2«.12*bi/ 


6«4tS«fif5 

9.5*«;7»7 


t  3i5,2.1,,8,4»0.,158.,85.,TU,!>8,,49i»34,,?5.'5»T9.,7<>,9<»0./ 

DATA  UPK7fl;60f6;72»5.75,5;42»5^0e»4.88*4.61«4.21*3.99«3iA7*3.8n« 


+  3, 7?, 3,86/,  DNW/492 . 0 ♦  285 , 5 ♦ 1 55 * 6 ♦  1 2T7t) , 9«> . C,T7v9 f 60 . 5 « 38 , 9 ,2:7. 7 , 
+  21.9*18,7,16.6415.1/,  CSW/30.7a9i0»13.8*i2i6,l0i6,9.69»8.35»6i48 


DATA  SHIELD  /20 . ,1 400. f?000,»3000, •5000. *10000; *20000. ,60000. / 
-  BATA  .  RATol.O  74. ,9,  ,25#  *49,  » 100*  e225»  «A00*  «625*  »6400./ 


Cffi)  *  ( •5+SlGN(<5*l«4Gn REAP (»An INI (20# tAMAXl (GTO# ) #»?) ) 

COMMENCE  .  . .  .  .  . ' 


READ  (L A 1 4 1 )  YIELD  '  '  '  '  . 

WRITE  (LBV41)  YIELD  *  . 


READ  ( LA y 4 1 )  HEAT-  ‘  ' 

writf  (lp»41)  heat  •  ...  .  ... 


READ  ( *; A 1 4 1 )  DIRTY  "  " '  '  ~  "  '  ” . 

WRITE  (18, 4U  DIRTY  .\  . 


'•  READ  MA,42)  JAIT  . . 

•WRITE  (IB»42)  JAIT  ...  : 


8 Y  r  YIC*»0.,333 
TY  =  Y!E«40>420 


DIRT  :  Oi  •  .  . 

,  IFUvEOiU  ;  DIRT  =  ,(2«9YiEf DIRTY (KW> )  *40*333 


DIRTP»i.5 

.1,18  .........  ...  .  • 


11  PSin(LtKW)  s  (BYOPSIRIUJ))** 


12  CALMUm  a  (TYttCALRCf.t - 

.  no  i  3  l  =  l  *  13  .  _ ..  ./ _  -  .  ■' . . 


J;  =  27-L  •  "  . . .  : 

RFNXdiKW*  a  »UPW(L)«nrRT  .  ...  •  . _  .  .  "■  •  -• 


RFNY(UKW)  a  «-CSW(L)#RAlN  ‘  ‘  '  • . 

PENX ( J»KW)  s  RENX il*KW)  ♦.  DNHtLI^RA! 


13  RENYU*KW)  a  *»RENY  (U  «KW) 

10  CONTINUE .  .  . 


20  REAn  (LA*43)  NAME » S TutNlTTrACTOR *  Y 

IF(  "EMS.lTiO.  J  GO  TO  24  .  . . _ _ 


0  ?r>  KT-  a  1,NF  . 

IF(.  NAHE.EO.NAMtYPUT)  I  GO  TO  24  . 


Nil 

GO  . TO  QO  . W  ' 


♦  44)  NAM£fSiOENS«r ACTOR *YTO*aTO*LAST ♦  Pal K  ’.fCALMKLfRtNMK 


mms 

■Hi 


TMU.CALMKURENMK 


* 

VSI'M 

no  21  j 

0. 

1  ♦  3 


7VK ( J)  =  OS 

NP  =  MINI  (DENS«OENS+a«  100* U 


1  I F (  riEMS.LT.O.  )  GO  TO  26 

IF(KT.GT « 1 j  NTS  a  M?ZiKT-l)  ♦  i 

_*  . 

•  '  *  , 

NTE  a  MF? (XT) 

WRITF  (LR*401)  NAMTYP(KT) 

% 

» 

PEV  a  -RAriT(KT) 

IF (RADT (KT) .LE.0»J  60  TO  23 

• 

C'x  Vs  * 

.  *_  _  ^ 

NP  s  MAXO(lfNP) 
AN  a  NP 


EV 
NPY 


N)  ♦  i 


icfciflfja 


u 


"WY  a  NPY-l 
WY  s  WY«0EV 


WX  a  WXttDEV 
.  CNEV  a  i.54WY 


FtnENS.GE.O.  V 
WR I  TP  (LAt47> 


YT  = 
XTO 


L  A5T 
I.R  a 


NTE  s  6 
YT  s.  YT  -  60 


GO  TO  23 


CrEV 

•CDEV 


:  XTO 

' 

-  1  « 

"  •  C  < 

! . CUS  C YT/RADIUS) 

1 

WX  =  WYREX  “ 

HFVX  =  DEV«EX 

i4 

REAn  aA*45).  STRING 
WRITE  (liR*45)  STRING 

• 

‘  -  , 

■ 

JR  s 
JR  s 


IR  a  LR  +  1 
WORTUl)  a  STRING (JR) 


rt*  >  **!*'<£ 


X  =  7-1  ■  ■ 

WORH<K)  =  ANn(WORO< V) .631 
LETTER (1)  =  LETTER  tn  /ftV 


I F  (  WORnTlT.LT  .0.)  WQRn(t)  S  "OR  <“WORO ill  I 

no  ?  r  s  ti6 _ 

I F  (  irFTTFRm  .E0»4fl  )  LFTTER ( I )  =  -1 
IF  (  LETTER  (!) . EQ.59  )  IFTTERU)  a  42 
IF<  I.ETT.E.R  <  I »  .EQ»60  VUmfM  n  a  10 
t.F'f  LETTER  ( I )  .GT.60  )  IFTTERU)  a  LE 
TFf  LETTER  (IK.GT.4f T“ UtTERU)  s  l£ 

tf(  letteru) ,ot.31  )  letter i  1  >  *  le 
TFT  UttEMiI)  .GT.15  )  LETTER U)  »  LE 

CONTINUE _ _ 

no  30  Tt  =  NTftiNTE  ' 

XT  =  XT  ♦  601  _ 

no  34  J  S  1*3  ~  ' 

I F (DENS«L  T « 0 » 0)  KLO  (I  T  •  J)  aLETTER  UT) 

AVEHOS TT)  a  0. 

PK (J)  ■  0.  _ 

IF (  S.GT.tM  READ  (4)  LBL (4)  • (PP(KI 

VALUE  s.LFTTERUT) _ . 

I F  <  npNS»i T«  0«  )  GO  TO  27 

VALUE  a  V A  Cl T ) _ . 

1F(RADT(KT).GT.0V)  GO  to  SZ'"’ 

NP  a  MlNUAREAUT)<»»l»»00»n 


-1 

42 

10 

LETTER  m*3 

rmrniH 

LETTER  U)«3 

LETTER  nr -3 


L  (4) • (Pp(K)  «PV(k)>iK«iV9) 


IF(  RA*T(KT) , Ef)»0.) 
NP.  s  MAXO(liNP) 
AN  *  NP  ‘ 

IF(  RAPT (XT) »Efl.Oi) 
NPY  =  SORT (AN)  ♦' 
NPX  «f  SORT  (  AN)  ♦ 


WX  s 


DEy 
•  99 


MINI  (DEN5«AREA UT)>il»400i  l) 


SORT (AREA (IT! 


WY  s 


WY  = 
WX  a 


chev 


WY^PEV 

WX«-OEV 


a  •  5<*WY 


EX  =  COStVUTJ/RAMUs) 


YT  = 
XT  = 


HFVX 


Y  (IT)  - 
X  (  T  T )  - 


=  DEV 


CDEV 

tDEV/EX 


IF (  VALUE. LE.-Oi 


0  33  TP  =  1 .NP 
DOSPSI (IP)  = 


GO  To  30 


4  . 

noscAuiP)  s  i 

33 

nnSRENliP)  S  i 

NWE  s 
DO  40  XW 


B  5 
E  = 


F(JALTflCW)  .LT.O)  fin  T 

DIST  s  SORT  ( AMAXl  (CALD  ( *P3  JD(i*KWL 


=  A*AXl(DISTf-*ENXm 
SOUT  a  AMAXl(DIST*-RENY(lfXWH . 


WEST  r  TAMAXl(DlST»RENX<26»KWr 
YNOPT  a  YT  ♦  WY  «•  SOUT  .  ...  ...  .  . 


XEAST  =  XT  •+  (WX+EASTiY'X 
Y SOI  IT  a  YT  -  SOUT  ' 


XWEST  s  XT  «  WEST/E 

no  60  !W  =  NWBtNWS  ...  .  .  . 


PWKtIW),. 


;  ■»  t,  t''  *  ’  -f 

c  ,  <■  v  -  ’  ,,  ^  *  t,  o  <*  *3  '■ 

A  *  2  M&ft  o'"  Z  ”  '^.'v  ’ 

0  o  .  v  &  v  xi  *'  a  Sfrte 


TF<YW<TW) .LT.YSOUT)  60  TO  60 


‘  IF(YW(IW)  *CT.-YN0RT)  GO  TO  60 
IF(XWMW)  .LTiXWEST)  00  TO  60 


I 


..  LX 


0  70  JX  =  1»NPX 
P.V  r  YT-YW(IW) 


LV  = 

if{lx.ot;26)  go  to  73 


F{PX.lT,RENX<lXfKW) )  GO  TO  72 
X  =  LX+1  . 


GO  T0  71  . 

LY  =  1., 


IF (LY»GT .26)  GO  TO  83 
IF(PY.LT.RENY(U*Xwn  Go  TO  62 


LY  =  LY+1 

GO  TO  8)  .  .  . 


LXY  a  MI  NO  (LX  4  20-LX  «L  Y«28**L  Y)  ' 
npSRFNMP)  s  DdSREN(TPJ  ♦  RENfLXY) 


PX*#2  ♦  PY*# 
ltlfl 


lF(niST,LT.CALD(L«XWn  GO  TO  84 

nosCALUP)  s  AMAXl  (DOSCAL  (IP)  «CAL  (L)PHEAT(KW) ) 

~  GO  to  85 

84  CONTINUE. 

!  ,''u 

35  Do  A 6  L  =“1,18 

' 

f  ' 

* 

.  IF (HIST.L  PSID  (L  «KW) ) .  GO  TO  66 

.  0 

/  ‘ 

,  i  <-  _  - 

86  CONTINUE 

.  87  .  IF ( IP ♦ FQi NP)  GO  TO  60 

*•  -  »  '  ‘  c  ’  0  : 

80  PY  s  PY  +  DEV 

,  70  PX  =  PX  ♦  DEVX 

> 

.. 

• 

60  continue 

50  CONTINUE 

no  '35  TP  =  l.fNP 

AVEOOS ( 1 )  -  AVEPOS(l)  ♦  tiOSPSKIP) 

'  '  ‘  11  r/ . 1  ‘ 

?  •  . 

AVEHOS  (?.) 
AVEP0S(3). 


AV£D05(2) 
AVEOOS (3) 


DOSCAt (IP) 
DOSRENfiPJ 


UPVAL  =  0. 
DO  51  X  s  1,9 


SURVAL  9  SURVAL 
PK  ( 1 )  s  PK(1)  ♦  ] 


PSFUN 


/RATTn(K))#PV(K)/VALUF 


~  "px  (?)  "a rpi(j2T  +  SURVAL  r”"“ 

SURVAL  S  sOrvALppSFUnC (n03CAL(iP)/eALMkLl 

pK  ( 2 F  =  PK  (2)  *  sIirVaL 

PK (3)  5  PK(3)  +  SUrVaL . 

~  .  .SilRVpF  9  SURVAL  -  T«  ^  ‘  •  .  “  r“" 

J  =  M T Ml  ( 39 • » OdlREN  ( I P )  /RENMKL. Kl  ...  .  M  . 

SURVAL-*  SURVALPPSFAlUJ)  ' 

in.  svLe.o.)  go  to  54 . . .  .......  . 

^  ‘  '  «  '  '  !j  '■ 

j  -  SURVAL  *  0,  •*  .  '  - . v  '  ■;  '  •  .  ■  •  '  - 

I  DO  53  X  s  jus  _  _  •  t . ■:  ■  . ,•  .  '  .  -  '• ,\  • 

! 

SURVPF  S  SURVPF  ♦  PF«K) /VALUE 

%  \  *  -  ' 

* 

-  -v  > 

•  V  *  '  • 

LI  5Uiffypr»LE«Q«K  go  Tn  53 


J  =  MINI  (39*  ♦  nOSRENMp)  /SHIELOIKI  )♦! 
SllPVAt  =  SimVAL >  SURvPFPPSFALLU). 


PK  U)  =  pk  13)  -  S 
39  CONTINUE 


0  ^  J  s  1  ♦  $ 

AV£nOS (J)  s  AVEhO$(Jl/AN 


PK  (J)  =  PK  (J)  /AN 

IF<  r>£MS.G£.0.  j  GO  TO  ,3  ft ,  * 


KLDHT'g)  -a  PK(1;)#VAlUE  +  i5 
KI.DUT»3)  =  PK(3)»VAlUF  ♦  .5 


3fl  ZVK(J)  =  ZVK { J)  ♦  pK(J)«VALUr 
_ VALUE  =  VALUE  *  (li"PK  M)«Pkl2)«PKIS)1 


VSU«  s  VSUM  VALUE 
IF C  DENS.GE.n,  )  GO  TO  39 . 


R L U ( I T • 1 )  s  VALUE  *  i 5 
GO  TO  30 


39  JT  =  IT-NTB+1 

WRITE  gF,4Q?)  GUAN (IT) , VALUE iPKiAVEbO 


IF(  LASTiEQ.O  )  VA ( IT)  «  VALUE 
30  CONTINUE  . 


IF «  PFNS.GE.OT)  GO  TO  91 
IR  =  LR  -  300 


DO  3  J 
IR  = 


28 


KLn  U  *  J)  *GTi*T4  ) 
IF (  KLP.M ♦  J)  iEOi.10  J 


klomtjf 

kLOltiJi 


♦  J)*4 


NUE 

MAP  f  IP)  s  KLliUtJ) 


(IR)  s  MAP  <  IR J 
IF (  KLD(ltJ) ,GT,3i  ) 


ONTINME 

IF  <  JRiLT,10  )  GO  TO 


IJlilifl'J] 


Hap  MR  j  .*  -MAP(IR) 


0  9 J  =  ii 

NT8  s  NTP  +  400  .  * 


s  NTp  +  LAST  -  1  . 

JJ  s  ZVKlJ)  p  FACTOR  . 


9  V/RilTF  (LB» A6)  J4TITLE  W)  i  (MAP ( IRI  ♦  IRsNT8»NTE) 

.  ,VSUM  =  V SUM' «  FACTOR-: — — - — “ ' 


WRITE  Cl)  LOL  U) #FACIL 
RETURN  .  .  ...  ,  . 


RMATfftH  CARDLBL * 12F6*2)  '  "  ‘  ' 

42  FORMAT I8H  CAR0L8L *1216)  . . 


43  FORMA  Til  XA6,A1«2F6 , 3  ♦  2F*  •  2  i-1.6  «3FJ1«1 


__^KF0RMAri(7lXA6«Ai.5X2Fi;3.5X?F6i2*!6«llX3F>ii:i 


"  45  FORMAT (4H  *10A65  .  . 

46  FORMAT  (14Hj  THIS  IS  MAP  «T1-«4H  OF  «A6// (30XiOA6)  j  -  ... 


47  FORMAT  (•/ j  8M  THIS  IS  INPUT  MAP".//) 

401  FOR.MAj <?H0  AA,SX5HVALUEl4Xi5HFRACT!0N  KIUED?IX1?HAVFRAGF  DOSE/3X4 


+HSITF7X9HSIIRVIVIMG1 1 X3HpKR4X3HPKT4X3HPKRl4X3HFS!flX3HCAL6X3HREN///) 
402  FORMAT  f?xA<5>,  ?XF7*  1»9X«3  j.3XF4*2)  i7X»3(2XF9«ll  «8XU) 


* I9H  ♦X! 


v£  >* 


SIBFTC  SHL«DK  DECK  • 
SUBROUTINE  SHELL 

common  ./PAR/  LA«LB 
CoMmon  /mOn/  Toen 
COMMON  /FAC/  NF,NA 

Wt-A  t /{TmjTnTwiwfT 

COMMON  /SHEL/  ZPFC 
HTRFN570R  PF (9) ,PV 


/MON/  I  DENT (fl) ,KEY ( 12) iLBL (4l - 
/FAC/  NF,NAMfYPC30),MFZC30>,VAl.Tt30)*RAnf«30)« 

rMirTT^mosTTv^rri^^ 

/SHEL/  ZPF(9,300)q?PVC9,300)  ,  _ 

ION  PF (9) »PV(9 ) ,5hL(9,300,2)  ,NCAt C2)7W0ff TUT 


QUAN (1000) 


EQUIVALENCE  <SHL,ZPF) 

“OAT A  NfAT v WORD ( 1 )  /8,9,l2H 
COMMENCE 

D0  sn  KT  si, NF 

IFj  IDPNT (5) «£Q,NAMTYP CkT) 

50  CoNTINUF  ~ 

51  NUM  s  MFZJKT) _ 

if.na.GT, i>  num  =  num-mfz ( 

_ in  MOD (KEY (2) ,2)  aEQ»Q  ) 

-  Do  52  J  S  l f NUM 

52  REAn  o)  LBl  (3)  ,  C7PFC*,J 


TO  51 


NUM-MFZ IKT-1) 

EQ  »0  )  60  TO  .60 


(7PFU»JUZPV(K«J) 


60 

GO  TO  70 

n0  41  J  =  NUM 

• 

DO  =  3,9  ' 

zpF(KtJ)  s  o;  .  ..  .  .  ; 

62 

61 

ZPV(K,J)  «  0. 

CONTINUE 

„ 

*•  ’a, 

1. 

63 

READ  (LA, 41).  KC»  J,  (PF  (10  ,Kal ,9)  ,NAME,LOC 

WRITE  (LP.,41)  KC»J, CPF(K)  ,Kal ,9) »NAME»LOC 

: 

64 

DO  64  K  s  1,9 

7PF(K,J)  s  ZPF (K « J)  ♦  PF (K) 

.  4 
|  4 

65 

IF(KC.LE.O)  GO  TO  63 

READ  (LA, 41)  KC,J,CPV(K),Ksi,9),NAME,LOC 

■ 

|  * 

WRITF  CLP, 41)  KC*J,(PV(K) «K*1 *9) , NAME, LOG 

00  66  K  s  1,9 

< 

■ 

I  A 

66 

ZPV (K ,  J)  =  ZPV (K, J)  ♦  PVCK) 

J F  (kC  »LF*  0.)  GO  TO  65 

_ =_ 

PFJK)  *  7PFCKi.ll 


TF(KFM8FO*0)  go  to  ta 
DO  73  K  5  KFM*fl 


s  5HCN  ♦  ZPFCKf.lV 
SHCN  =  AMAxi,(FRACFPZPF<9*  Jj^SHCNtOil 


D 

88 

0  88  K  s  2*N 

4  s  A  ♦  SHL (K» J«M) 

", 

- 

- 

A  s  AMIN1U,,V#UTIL/A) 

00  89.  K  s  2,N  _ . 

R 


81  CONTINUE 

90  IF(  KEY (2) »LT»2  ). RETURN 


LBL (2) 

no  9i  j  s. 

RETURN 

41  FORMAT ( 1H  »I1»I4*9F7;HA6.A4)  ...  .  •  .  . 


FORMAT (8H  CAROL 8t 1F6# 3»? (AX fF6#3» 16 « j 6) ) 

44  FORMAT (40H1  FALLOUT  SHELTERS  BY  AREA  ♦  CAT  9  *  POP///IOH.  CATEGORY. 


wWffwft  iFFlf  if  [ft tf1!  i  \  4At  l 8  1 rr« k  1 4 1  T/j*m 

■M/ilim  f  -01  V.YJr.T;U.WS/J  TiTTWC*  \  m 


+  5XIHl9XlH29XlH3qXlH49*lH39XlH69XlHT9XiH*9KlH9  //I 
7  F0RMAT(?4Hi  sHETTER  SUMMARY  TOTAL  $T?5X»4HPOP  FlO# 1///10H  CATEGORT 
♦  5XiHl<?XlHR9XlH39XlH49XlH»>9XlH69XlH79XiH89XlH99XnHc6NSTRUCTED//:. 


♦  10H  FALLOUT  nr  IO«tl  *  l“Xr 10*  j/FlUn  0LA5T  ?no#H _ _ _ _ 

+//38H  CONSTRUCTS  TO  BE  ADnEO  TO  CATEOORf  EX«  213*2X1 2HftgSPECTTVEllV) 


FORMAT (30H1  SHELTER  UTILIZATION  tPARAn  SF6*3« 

5  X 1 H 1 9  X 1 H2  9  X 1 H  39  X 1H49  X  t  HS9  X 1 H69X IHT9  X IHBIX IHT// 


*A- 


& 


u_. _ _ 


.  _ _ _ 


i  j*«~ 

*  Ji— 
n~~ 

v _ 


SUBROUTINE  SURVEY 

'  0 1  HENS  I  ON  S  V  2 1  •  0  U2 1  •  Old  I  ) ,  Q3! 1 II  ,  05 1 1 1  IYP  111  I*  BACHN I  nT.WJ25)T* 

♦  T8STG( 25 )  ,Tc$TGI25)  »THRH25),THR2(25I  »AZlTH(25l,ELETHf25l  *NNC25>  r 

♦  ANOZ  <  2  5 ) , R AtM ( 25 ) ,COI 25 ) ,0C  til , 25 1  ,WT<23)  ,TCR0S5<3 I 
DIMENSION  RADARt  8) «XSf  8 1  o  YS( 8 1 • ZS( II • 

♦  Dl  l  { 8 )  ,012  ( 8  >  ,  021  ( 8  > ,  0221 8 )  ,02318 ) ,  Ml  HI  ,0921  81 ,033  m - - 


o  O-  ' 


r 

•«— ■ 


-  RADF(X)  *  RA1  ♦  RA2*C0S  (X+X)  ♦  RA3*C0S  UXmMX+m: 

GEOFtX)  *  X  ♦  F1*SIN«X*X)  ♦  F2*SIN(t X+X J  +  1  X-fXM 
CONSTANTS 

DATA  BACHN/0.0,0.8, 1 .0, 1*2, 1*5, 2*0, 3*0,4, 0,6*0, 6*0, 10*  0/ 

- -  DATA  S  /11.63, 2*9. 5595, 2*10.9169,7*8.3585/ 

DATA  Q/-5., li. ,25. ,47. *52. , 79. ,90. ,10 5. , 160. ,170. *200. ,700./ 
qATA  Ql  /“  .65807881, -1.0105494,-3. 2178026, -6.5U0337»-7.2373653,  ~ 
♦rlO. 759466, -13.028863, -15o977971,-20.7732T7,-2l.094464,-21. 878067/ 

- -  DATA"  Q3  /1-.0Z02628,9,-'* 0000131 8m013846579*V0000Z35B4Y-.  01 593 154, 

♦.000032939, ,024145841, ,088628908, .007543412, ,003507148, ,002221291/ 

- —DATA  05  /-4.2581444>-n95ei974V"12.383264,'3126.0905T=6.5856931,““ 

4,6264.4021,  9.5411985,  2.7082397,  4,4164793,  7,8329587,  10.761370/ 

''™V"  DATA  P/3709. 631, 474. 088, 51^7632, 2.64295,  lo29325,.026Tf4V 
+.28265E-2,  .2049698E-3,  .839132E-5,  .6585E-5, .  3402E-5/  ._ 

OATA  PI.OEGREE, OMEGA, AG/3. 1415927*57. 29578, .0043752695, 19.061428/, 
♦.  G,AE,Aj,AK,H  /22SSQS040.  .3443.939,  .32469E-2.  .54E-4,.1075E«?4/, 


U-., 


+  RA1,RA2,RA3,C8  / 343 8. 1533,5. 7978506 ,-.i 2222 54E-T,23. 277^ 

♦  R1  ,R2, PHO, SPO  /6355. 766,3432. 3791, 21 16. 2, 1.48121/ , 

*■  H I  NO ,  BELOW ,  HFVH,  H IGHT  V-1234.Y-3*  ,“29.~624*~52.7 » 

♦  L A  j  LB , LC  , JPR I NT , T INT ,  TPRT , T$T OP  /9, 6, 4,3,. 01, ,1 ,100,/, _ 

4  LOCI ,  IC,  IPAGE* 

COMMENCE-.  . 1  .  ..  _ _1 - __ - - - - 

F1  -  -.337276-2 
F2  «  .567E-5 
REWIND; LC 

RADIAN' «  l.  /  OEGREE 
lOo  REAO  (LA, 45)  ’ 

LAUNCH  *•  0 

9i  i page  «  i Page  4  1 

LINE  *  0  ■ 

WRITE  (LB, 45)  IPAGE 
CARD  INPUTS 

1  READ  . (LA, 48)  K,A*B,C,D,J 
IF(  .K.EQ.O  )  RETURN 


'J 


GO  TO  1 80% X *  82, 83,84,85,86,87, 88* IfltK 


•  ,  »n*' 


*  * 

* 


•  'f 


2-27 


•WP*': 


)  >  t. 


<t 

.  -  Cf 


>  *  ^ 


Y  ■  '  ^ '  S’&fie  c  -/■  j| 


t^fcWMwiWiii  m  imwm  wiiwiiiii  *■»»»'■*!>*  «>  —  rwriwr^ 


86  CONTINUE  ''  V*"  .  ••"'"o’™  ,  ■  ‘  7.  I* 

..  .  .  89.  RETURN  .'  J  _ 

CARO  0  BROUGHT,  IN  RADAR  LOCATlONjl  .  .  ;  :.  /‘ 

•  '  8C  IF!  J.E0.9  )  LOCI  «  6.'/'  ’./'V  -  u  .  ;  .-..  .  - 

t-OC.I  *  LOCI  ♦  1'“;.  rT ;  '  ** 

■  GLAt  ?  A  *  RADIAN  ,  ••.•:•  _ - - - - — ....— : 

FLON  *  :B  •  RADIAN  '" .  ... 

UP  *  COS  (GLAT)  ......... _ 1.^^,  _ -w-J— . 

•  '*  UZ  -  S IN i GLAT)  '  '  ,  , 

CFLON  ■«  COSIFLON)  _ „ _ _ _ _ _ _ _ 

SFLON  ■*  SIN(FLON) 

.  •  FLAT  *  GEGF(GLAT)  _ _ *  *»«*,*.. _ — - - 1 — . 

RE  -  hADF(GLAT)  :  ... 

HOR  «  RE»C0S3FLAT)  »  C»UP . _ _ 

CARTESIAN  COORDINATES  OF  SITES 

ZS(LOCI)  -  RE»S1N(FLATI  ♦  fUZ  .  .  . 

xsaocil  »  H0R*CFLQN  .  ...  •' . . 

YS(  LOCI )  «  HOR»SFLON  .  ‘  ...  . _  *_i.  • 

;  RADAR ( LOCI  )  ”f  0 .  ~  T  '  ~~  Tv? 

COORDS  EAST#  NORTH  #JJP  AT  SITE _ _ _ u _ _ 

01 11  LOCI)  *  -SFLON  ;  1  \  i  ■  ■ 

012<L0CI)  «-  CFLON  .  a..':  •  '  - ■  .  •  .  ■  ‘  '-  •  : 

021  (LOCI )  *»UZ*CFL0N  . 

D22(L0CI)  «  -UZ*SF10N_  .  ...  ,r.  .  .. 

‘  023 (LOCI)  VrUP . I*  ““  “  *“*'•  '  '•  *  •  ~  ““  ’  •7““  “ 

D3KL0CI)  *  UP*CFL0N  .  _  J 

. .032 (LOCI')'  -  -WW» W - -.----r' - . - ; - ?»“ Tj 

D33( LOCI )  -  UZ  '  . 

- WRITE  (LB.AZTTOCTilJBTCinr  “  — — — :  ~ 

CO  TO  1  ‘  -1  • _ .  •  \ 

81  "  JSAT  *  J  -r-*—-- f- 

TINT  «  A 

"““TPRT-B  ‘  ~~”7 

\  TSTOP  «  C  . -  . 

- - N0RBIT~*~D~  - rTr-"“  — “ 

I  BIN  *  1  -  JSAt  .  '  _  ;.L., _ _  IL-;,L  ■  .<-.••  • _ _ 

GO  fcTi  •  -  “““■  “  ‘ 

82  IF(  JPRINT.GT.3  )  GO  TO  21  '  -  -Mnwi- 

CARO  3  BROUGHT  IN  INI' Ml  AND  PINAlVAlW  WTHRUST  AMD  YlNf  TdR  ^U8SttGn 

.  83  IFI  J.E0.9  I  IC  *0 _ _ -  . '-.....-V: _ : _ _ _ 

ic  *  ic  ♦  i  *• 

THRU, IC)  *  A  .i.7.  .  _ _ _ _ _ _ _ 

THR2MC)  ■  B  .  ’ 

TBSTG(IC)  -  C  _ .1 _ _ _ _ _ _ 

TESTG(IC)  0 _  ' _ •••  ‘  '■  _ _ _ 

WRITE  1t5t9t0nCiwtBt£fW  ■  r(  -  -  --  ••  - - - 

‘  '  ‘  GO  TO  1  .  . 


♦  *  . 


>** 


A  ,  «! 
»  J 
*  S 


CARO  A  BROUGHT  IN  SPECIFIC  IMPULSE  AND  THRUST  OftlEN*  VTION  FOR  EACH  SUBS74GE 

84  NN(IC)  =  j  - - — - - - - — ~ 

RATH(IC)  *  A  *J. 

A2ITHUCI  «  B«RA01AM  - - - - - 

ELETHdC)  -  C*RADIAN.  ; 

ANGZdC)  «  0  - - - - - - - - — 

WRITE  (LB. 911)  A»0,J*C,B  ..  . 


CARO  6  ENTERS  CO,  AREA  ANO  WEIGHT  FOR  EACH  SUBSTAGE 

85  CD(IC)  «  A  - - 

_88(IC)  *  B 

WT(IC)  *  C  ^ - 

.WRITE  (LB, 912)  A,B,C*J 

IF(  CO(  IO.GE.O.)  GO  TO  850  --w 


_ -CAROS  5  ENTER  DRAG  COEFFICIENT  TABLE 

RE  AO  ( LA , 43 )  ( DC  *  I.IC  )Vl-lt"lf)  - - - 

WRITE  (L8.913)  BACHN. (0C( I. IC>V I -1.11) 
850  1 F (  J.EQ.O  )  GO  TO"l  '  V' 

ICC  -«  IC  +  J  -  l 

00  330  JJ-IC.ICC  “  ■  — 

C0(jj4l)  «  CO(JJ) 

BB( JJ+1 J«BB( JJ )  “  “  - 

IF(  CDdCJ.GE.O.)  GO  TO  320 

00  852  i'm  ltil  . . — ~ 

DC ( I , JJ+1 }  -  DC(I.JJ) 

'330  CONTINUE  . ~ . - — — - 

GO  TO  1 


V - . 


852 


•o  ! 


c  , 

•>  N 


card  i  brought  initial  LATirabrrroNsirToerrwciemTkiiirTiiir' 

87  T*D  >*'  '  4  •'  *  • 


r  ------ - 


IF(  J.NE.2  )  GO  TO  8t0 “ 
WRITE  (LBt601)AfBtC,0 

X  »  A  *  - -------- 

Y  *  B 

Z  «  c  - - — — 

RSQ  »  X*X  4  Y*Y  *  2* Z 

R  »  SORT  (RSO)  . . . — 

HOR  «  SORT  (X*X4Y*Y) 

FLON  *  2*«ATAN( Y/CH0R4JO )” - ~“ 

GLAT  «  AT.AN  ((1.  4  C8/R|*2/H0R> 
GO  TO  872  ^ 


•)MN 


- - ( 


870  WRITE  (LB, $01 )  A.BtCfOfJ 

■“  '  GLAT  "  A  *  RADIAN  - - . . 

FLON  *  B  *  RADIAN 

i f<  j.le.o  )  FLorr r*t»tJMEGir" 
872  FLAT  »  GEOF(GLAT) 

RE  «  RADF(GLAT)  - - - - - 

UP  •  COS (GLAT) 

- UZ  »  SIN(GLAT) - : — - — • - : - 


W..—  —  -■  -  —  r-^T,  1 - — 

t  ,  '  "I" 


Mil. _ MAJI  J  Mi  All  A, 

wrtun  gusaruuna 


V, 


,  o  - 


SFLON  «  SIN (FLON) 
IF(  J.EQ.2  5  GO  TO  1 
HOR  *  RE«CQ$( FLAT ) 
Z  »  RE*SIN(FLATI  4 
X  -  H0R*CFL0N 
Y  «  KORoSFLON 


4”  C*UP 
C«UZ 


— — — — — -T~— — — — u— ; 


2-29 


...  -  s>! 


&  J- 

J  'o  ”  ?} 

</>  '  _v  *  <*<>'<>  \ 

v  •  "i  .  <  Q  * 

s<»  -  4  ;  ■=>  „  •  '  . 
c-  v  .  ^  •  Of  * 


IL 


•  ^  O  v  '' 

'  V  ,  V  *  ,  2 


■•AM 


.  ,  0  ■' 


-.t  .  ^m:{. 


4 


GO  ,iu  l  ■  ,  . 

CARO  8  BROUGHT  INITIAL  VElOCltY  IHEADING*.  INCLINE?  ANO  tfAGNIfUDEJ  _ 

88  1F(  J«NE*2  I  GO  TO  830  ~  " 

WRITE  iLB»602|.A,BtC  ■  ...  -  ...  ;  -  v.  .  .  -  . 

VX  =  A  '  ’  "■  ‘  I  ~  ~ 

vy  =  b  •: '  ..  •  v  -  .  ;  .  . 

VZ  a  C  ~~  .~rr~T  .  ;  -  77- 

GO  to  102  •  •■  .  ..I  ‘  .  I 

880  WRI TE  ( LB  ,  902 )  A  ,  8  »C  ,  J 

A*A*,RAD1AN  '  .  „  .  .  . 

' 8»8»RADIAN  . 

U3  »  COS ( B } *C 

U1  «  U3*SIN(A)  “  ' ------ ' **—“ 7  “ 

ti2'  *  u3*cos(.Ai  ...  *  •  '*  •  ;  • 

U3>  «  SIN i'BI«C  ■  ““  ~**r~* --e— 

_  VZ  -  UZ*U3  +  UP»U2  -  1  .  '  ■  .  .  '  .  .  ’  :  . 

U2  «  UP*U3  r  UZVUZ  ““  *  ; 1  “  *  :  .  “  7  7 

VX  .a  U2*CFL0N  -  Ul»SFLON  ..  ,  .!■ _ 

VY  »  U2«SFL0N  ♦  Ul»CFLOM  ,r  ^  J  >  ~  ‘  •  " 

IFI  J.6Q»1  I  GO  T0  40.2  '  . 

VX  -  VX  -  Y  *  OMEGA  ..  ~  ~  “  ~  ^TT?  --r-r-^-rnT' 

•VY  *  VY  ♦  X  •  OMEGA 

C  BEGIN  LAUNCH  ~ .  “  .  .  “ 

..  102  CALL  SLITEi II  _  _ '  J:  ‘  ... ,  : ' .. . 

LINE  a  60  T  Tr"T“  :  ‘  •““*  r  ~ 

LAUNCH  «  LAUNCH  ♦  1_  *  -v  ‘ 

LEFT  «  0  .  ‘  ~ 

1  l i st.  *  0  .  . 

IH  *  ii  ; . ‘  „  "  *  ■  ••  /  :  r“  “  ^ 


ASSIGN-  705  TO  MM  ‘  ..  _ _ _ _ 

l  *  1  .  r  rr  ■  • 


NNN  «'  NNILI 


W  -  WTdi  •  .  ,  •...  -777^-——,—* - : 

tl  *  t  .  •  .•  .  *■  *  Ui 

"  ** TB  »  T  “T-----;--  -■  .  •  7“  !r  : 

TINY  »  TINT/10.  -  r  “  -r-r-T  ■  | 

CALCULATIONS  FOR  HEIGHT  AND  ACCELERATION  "  1  ,  , 

'  68.-  RSQ  «  X*X  ♦  Y*Y  <•  Ul  vr;_  - - - — -r-— ^ 

R  *>  SORT  ( RSQ)  ' 

'  SINO  «  Z/R  '  -  . .  . . . . “n  * 

SINSQO  «  Z*Z  /  RSQ  .1  » 

■'*  "  '  A  TO  ft  SQ  «  AE*AE/RSQ  . .  . . .  . 7'  *  4 

GR3  *  G/(R*RSQ)  ,'rf  “i 

COEFG  *  -GR3»AT0RSQ*H*U6-3.*SiHSQ8T»AE . .  ‘  '  —li  1 

COEFZ  ■  -GR3#2.«AT0RSQ»(AJ  ♦  AK»AT0RSQ*M.-l!NSQ0#V,/3. } | 

‘  COEFR  *  -GR3*L1.  ❖  A At 0RSQt4|T»r«VnaMS(»T'^AlKf 08^  f  | 

♦  SINSQ0*I0.5#SINSQ0«*2I  -  V~~ 

;hor  *  SORT  (X*X*Y«YJ  *'"T  “  - — - - - - - •>»  ^ 

FLON  »  2.«ATANIY/CHQRtXH  . 

.  cflon  » cosiftLo^ir  — -r— — ~r — 


>ro 


A  : 

■  4*1 

-jV  oS 


SFLON  *  SSN(FLON) 

“GUAT  "»7ATAN‘ -rriV’  ^CE/.RTpZ/KORT" 


.  UP  «  COS  (GL  AT  I  ....  ■; ^  W 

......  ...  .  .2^30...... 4.s ' .  ,  .  ».  .  .  -4, 


,  ^  ^  I 


a:  ^ 


if  V°,A/ 

ir  e-  "  vv""  *•.  '  ■ 


'Jl 


&  ,  , 

i  * 

W 


-  6 


'  UX  *  UP  «  CFLON 
UY  =  UP  *  SFLON 
UZ  =  SIN(GLAT) 
FLAT  ?  GEDFtGLAT) 
RE  =  RADF(GLAT) 


sjsfessssss§sg^afS^SS^S^i 


HIGH  =  UZ*(Z-RE»S  IN  if  urn  _ _ _ _ _ _ .,  - 


<at 


<*> 


XDOT  *  VX  <■  Y*OMEGA 
YOOT  =  VY  -  X»OMEGA 

SPEED  »  SQRt  ( XDOT *X DOT  ♦  Yb0T«Y0OT~«'  VjZWfl  ~'' 
U3  «  UX»XDQT  ♦  UY*YOOT  *  UZ®VZ 
U1  =  UX^YDOT  -UVwXDOt  ~  " 

U2  =  VZ  •?  UZ»U3 

U3  ®  U3/SPEEO  '  '7  '  “ 

HEADIN  *  2*  *AT AN ( Ul/ ( U2+SQRT ( UI*«2*U2**2 1 1 ) 


,  \ 


IJ4*  —  •**  -  - 


53 

54 


5t 


ATTACK  *  2,*  At  AN  (U3/  «1 .  VSQRf^T*^W*l»2T)T 
XF(  HIGH.LE.O.  J  GO  TO  71  ...  ■ 

ATCH  =  R  X  *H  I GH/tR2+fT  I GH I 

IFi  ATCH.LE.QU2)  »  GO  TO  54  '  *!•  \  .  ■ _ ft 

PH  =  0.  '  " 

COEFV  =  6.  .v  .  . 

■  GO  TO  70  T  TT-  “'“7"  “ 

IH  Si  TH  -  1  . _ 

IFi  ATCH.LE.QUH)  »  GO  TO  53  . 

IFi  ATCH. LE.Qi  IH+1 J  I  GO  TO  57  „ 

IH  *  IH  +  l  '"'7  “““ 

GO  TO  54  7  . 

PH  -  A L 0 G  (1 .  +Q3t  I H I  *( ATCH-Qd  Wm  “  1  “ 

DENSE  *  EXP(()1(IH)~Q5CIH)«PH)  _ ^ _ .  .. 

ph  -*  p  (  i  h  )  *  exp  i  pH«  <  r.-05  <  iH  n »  -  ■  T  •••  ' 

SPSOUN  «  SUH)  +  (SiIH*n-SiIHn«iATCH-0(IHn/(ni3H+I)-QiIHn 

BACH  *  SP EEO/ SPSOUN  . . . —  — 

CP  AG  *  CO(L)  . i.  ‘| 

IFi  DRAG.GE.O. )  GO' TO '59 .  - -*rr-r“ 

00  56  K  *  1,11 


‘  56 

58 

59 

”*  t 

736 

•  70 


IFi  BACH.LC.BACHN(K)  )  GO  TO  58  ' 

DRAG  *  -DRAG 
GO  TO  59 

BACH* { BACH-BACHN ( K-l  H  / i BACHN <  K  >“BACHN«  K-H i 
OR  AG  *.  DC(K-1»LI  +  <0C i  KYI) 

DRAG  *  DRAG*BB(l)/W  • _ 

COEFV  «  -1 . 31708i6*DftAG*DSNSr*5FEED“ - “ 

GO  TO  70 

SUBSTAGE  L  IS  OVER~'~~ . — — r-— 

ASS IGN  705  TO  MM  ‘ 

W  *  W  -  ;T  E S T G  i  L? -T B  *TL  )«THRUST/fCATH at - 

L  *  L+l  -  . 

IFi  WTILJ.GT.O.  )  H  »  TiT(L) - - - li““ - 


r 


— !i 


705 

710 

731 

-7<i7 


NNN  *  NN(L) 

AT  HR  «d.  •  . 

iFH.Gt.ici  go  to  701 _ ^  _ _ ;  •,/  .  ; 

IR  |-  I  L  .LT.T55TG  i  L  ilr  GOTOTvl-  “  _w‘  ’ 

IF(T-TL.GTi-rESTG(in  GO  TO  736  •  .  ..  7. 

GO  TO  MM, (705.732)  *  1  r“  "  ~ 

CHOICE  OF  ORIENTATION  AT  BEGINNING  OF  SUBSTAGE 

ASSIGN  732  TO  MM  - - - : - 

GO  TO (7 10, 712, 708), NNN  .  *  ■  •?*.•  ‘4  7.  • 

I F  (  SPEED .  LEiSPO  ") 

IFi  HIGH.GT.HFTH  )  QO  TO  708  •  '  ; 

THX  «  'X00T/SPEE0~7~^~7 - , 

THY  *  YOOT/_SP_EEO_  V.,  •  ;  . 

”7^ 

*  X  V_.  -H  .*  -  ^  »  *  *  %  ^4*.' 


■  ...  m»+mm  •rnmmmmm  ^  if  ,  In 


777-'  y 

.'V'X'  “v  '■?,  ■  v  77v  '*9 

^  .r  v  C  -  ;  *  V’  g  w  ^  7. 

'  -  *-  f/  <'*>v  . <*  ' r  ^7?  'S.  v  4  .k  ,,  , 

A  -  y.% 

■■■  -  ,» f  c  r*.;  •;  o  ■ 


m 


<  >  „ 

’  '  «  .  ;  ?V  '  %  ^  (y  -  ,  *  .  <* 

—  „ .  .  .  .  ,  _  - —  .  mMmmmIm  *•**£>♦»«**•  4s#  tiiMM  m*  «■«««•• 

THZ  *  VZ/SFEEO  ,  ,  /  r:'-' 

GO  TO  708  ‘  >  ■ 

712  113  *  COS(ELETHM))  *“•'  \  r~  ~~~ 

U1  i  U3  *  SINtHEAbiN+AZITHID)  ’  ..  IV  '»  ’".  '  ■  22  IZlH  v'--  -•••.■' 

U2  «  U3  *  castHEAdiN+AzifHfCn  '  "  ■’  -  '.  .  :••  'V",  • : 

U3  *  SIN( ELETH{ L) )  V  ;: 

THX  =  U3*UX  -  U2»UZ»CF10N  -'UI>SFL0N  . . ' 

THY  =  U3*UY  -  U2*U2®SFLON  ♦  Ul»CFL0N  .  .  ’ 

THZ  *  U3*UZ  +  U2*UP  '~7 .  .  .  *•"  “V" 

GO  to  708  _  ■ _ * _ •  ;■■•■.  ,  ;  7  r 

C  ORIENTATION  DURIN*  SUBSTAGE .  ~  ~  'V 

732  GO  TQ(733, 708,708) ,NNN  '  _  J  _ 

‘  733  IF(  SPE60.GT.SP0  )  GO  TO  731“’“  .  ““  . ““ 

703  THRUST  ■*  THRl(L)  +  UTHR2M)  -  THRl I LH / tTESTGC t I  -  TBSTGCLH  14(7-  ^ 

*  7L~7B5TG(D)  . . "7 

IF1THRUST.LE.O.)  GO  fO  701  _  _  _  '  '  r'' A.SS 

AT  HR  «  ( I.+  ( PHO-PH )  *ANOZ  C  Lf/AHAX 1ITHRML  )*“fHR2f(  L I  O'®  AG*  THRU$T/W 
W  »  H  -  THRUST*  AM  INI  (  TINT « T-f  L-TBSTGl  LS  ) /RATM C LI  „  « 

701  AX  «  ATHR*THX  .♦  CQEFV*X00f  ♦  X#C0EFR  ~ 

AY  *  ATHR*THY  +  CQEFV*YDOT  ♦  Y*C0EFR  _  _  j 

A Z  *  ATHR*THZ  ♦  COEFV*  VZ  "  ♦  T#“(CQEyR^TOEFZT“+  CQEFG  7“ 

CHECKS  FOR  PRINTING  AND  STOPPING  _ ..  .  .  '  ‘  ■  ■ . 

GLAT  «  GLAT*DEGREE  ...  .  f  • '->  y  v  | 

GLON  «  FLON  -  T  •  OMEGA.  _  _  .  .J.  ’  ..  ..  V  /  j 

XF(GLON.LT.-PI)  GLON  *  PI  ♦  AWdTGLW,~PIT7“7  ““““T  ^ 

GLOW  «  GLON*DEGREE  _  '  ■  V  ‘  ’•  : 

HE  AO  IN'  •  HEADIN*0EGR6E  •  “  “  ““""7“  7r'7~  >  '  1  •  :  i 

ATTACK  -  ATTACK*0E6REE  .  > .  .  .  ,  ,  .'j.......  ,  -1 

I F  ( T+T I  NY  •  LT  •  TREC }  GO  TO  64  "  .  ,  V  7  •  ,j 

C  WRITE  BINARY  TAPE  ONLY  ,  7.?  ..  } 

I F I  IBIN  )  62,60,61  ?  .^.-Trr-ww-TO-r-rl 

62  WRITE  aC»  T ,  GLAT ,  GLQN  ,H  I NG  ,X,Y  »  Z _ }  777,..-^>  !.  7  .~.f 

Jn^n^n  ’  *  *  ”  ^  7  ^  *->V  'V '  1 

GO  TO  60  .  • _ ,  ,i; 

“  “61  WRITE  (LCi  T ,  GLAT  t  GLON » HIGHVX*Y#Z  •  *  “  “  ~  .{ 

60  IFC  M0D(JPRINT,2I.EQ,0»  GO  TO  67  ,  .  J 

C  PRINT  TRAJECTORY  DATA  ’  ”““"“1 

ATOT  -  SQRT(AX»AX4AY*AYtAZ*AZ»/AG_  .  '  M 

ATHRG  «  ATHR/AG  _ 

DRAG  «  -COEFV ‘SPEED/ AG  *"7“  ^ 

.  IF1LINE.LT. 54»  GO  TO  66  “  ' ; 

I  PAGE  *  I  PAGE  ♦  1  •  .  ..  V  J 

LINE  »  0  -“7“  “  ~~  . “rr”":; “““I 

WRITE  (LB, 45)  IPAGE  , 

WRITE  (LB, 44)  —  -r— - r  T----TrTj--.r-«-H 

66  WRITE  ( LB , 47 )  T , GLAT , GLON i.H IGH* HEAOINtATTACXc SPEED* ATOT t ATHRG* DRAG  . , 

LINE  LINE'  +  1  “  '  '  •  "  >  '  >7-  “  “  •  ■'  <r  | 

67  IF(  JPRINT.GT.l  )  GO  TO  29  ,  .•  •;  ;  .  . '  .»  •  7.  | 

'230  TREC  »  T  +"  TPRT  ""  . - 7^ - ---—^-------,---4, 

64  IF(HIGH.CT.HIGHT)  GO  TO  19  “  '  7  Z*  ,  7\  :'4 

'  ‘  I F  { LEFT .  EQ .  6 )  ISO  TO  '206  '  - 

IFC  JSAT.EQ.l  i  GO  TO  71  '  '  *  ■  .7  .  '  '  '  '  ,  J  il 


ir \  ruvvn.ur.v/.  i  w  iu  .cud 
C  THE  END' '  \  }■' 

71  WRITE  ( LB, 49)  "L AU NC H,G L ATVS LONVTT” 
GO  TO  1 

19  left  «  i  - — — 7*““ 

206  IFC  T.GE.TStOP  )  GO  TO  90 

C  INTEGRATION — : - - 7 — - - ~~ 

TB  «  T 


'SSSSSSSSSSS 


■caUrii 


kfiii 


iiiwi>nan»w 


'.  s  .  •  '•  -V\-  "' 

•  •  •*  ?  *  v.-X-  '*  ':r  *  *  <k'm‘  * ^ 


**  *  ^  ■  A‘« 


"*‘‘j 


28  2 

T  =  VIND<T,TINT> 

X  =  OPNV  (  X,  VXJ 
V  =  DPNV  {  Y,  VY) 

2  =  OPNV  (  2,  VZ)‘ 

VX«  OPNV  <  VX»AX) 

VY»  OPNV  (  VY,AY) 

VZ=  OPNV  <  VZyAZ) 

IF{  JSAT.EQ.O  )  *  GO  TO  68  “““ — 4~- - - 

IF(  Z*Z8.Sf.O.  j  GO  TO  68  ’  *  '  >  *  ■  '•  /■«/-' 

Equator  crossing - *  — - ; — - — r— — — — —— ~ • 


i< 4  ’  ■  \\  •  ;  "* »  *  i  's  v‘  *  -*  */’•  1 

O.  * 


y  f 


Km 


J’  * 


1 


K- 


ILIST*IL!ST+i,  . 

I8IN=1{  IUST)  “  T  -  CV2/A2 |VaVrl5WY^ 

I F  <  I ST-3 )  6S  ,10X2  * 2012 - - - - - — ~~ — 1 

_  1012  IBIN=-1 

GOT068  - -  *  '  ~  - - -■—■■.  ,• - - - - - 

2012  IBIN=0 

GO  TO  68  . . . — - - - - 

C  OVERTIME 

IF<  INDEX. GT. 2  )  GO  TO  1014 - — 

WRITE  ( LB, 1015 ) 

G0T0499  . . _ 

1014  PERIOD  «  ( TCRQSS (INDEX ) -TCROSSi l ) )/( •  !•( FLOAT  ( INDEX 1-1 .01 1  ! 

WRITE  <LC)  PERIOD"  '■  ■  •  ^ _ 

004000 jUST*i,IUST  !  .  % 

4000  WRITE  ( L8 ,  4001 )  JC~| Sf VTCR61GrjQGT>^‘ ~ — ‘ 

..  WRITE  <18,201611  PER t GO  *'  '  ‘  •  . 

4.99  WRITE'  ( LB , 606 ) T - ■ - - 

go  to.  i  *\  '  '  ’  ■  .  ;K*\  \  cv/4 

C  TAPE  INPUTS  — - 

_ 21  READ  <LC)  T, GUT, GLON, HIGH, X*V,Z  ’  "  -  />  r  4 

T  =  T  +  SHIFT  - - - “ 


A  ‘i 


J 


— u_ 


©  { 


_  eARTH  (ROTJATION  to  present  time  ’ 

29  ROT=T*OM£GA  — . - 

_  COSOMaCOS  (RQTJ  •  ” 

",  "  SINQM=SIN  ( ROTil' - - - 

_ XP=X*CQSOM  +  Y»SINOM 

YP=-x*SINOM+"Y*COSaH — -  ••  ■  .  ..  „ _ T 

c  TRANSLATE  *■  '•*'[' 

~  DO  25  I  *  1,L0CI . — — 

xv*xp-xs<=n  '  1  ' 

YV=*YP~YS<  I )  . 

_  2v*z  ~zs<n 

PZ  a  xv*03i<  I )  +YV*D32{  II 4-ZV#tJ33'(TT 
IF {P2.LJ. -BELOW)  GO  TO  25 

C  HEIGHT  ABOVE  RADAR  PLANE  - - 

PX  »  XV*  DIHI)+  YV»0i2m 
VER2  *  PX*PX  +  PY*PV  '  - - - - - 

PY  «  XV*D2i<n+VV*022U)+2V»023III 
VER  *  SORT  (VER2)  - 

RH02  «  VER2  +  PZ*PZ  , 

C  •  RADAR  SPHERICAL'  COORDINATES  r~- ----- --- 

RHO  *  SORT  (RH02)  •  .  V  '  I 


<<  ' 


ELEV  »-  ATAN<Pr/VERT~#"OISRJE“ 

AZI  *  2 • *AT AM<  PX/ <  VER^PY ) )  t  DEGREE 

I F  <  L I NE .  LT 1 54  PGO-TO-JA - ‘~s- — 

I PAGE  *  I PAGE  ♦  1 


— . •  '  .  ;  .  (■ 

. . 

. .  ,  , ,  .  .  ? 

7-vri; 


4  ■i  »  . 


.SLl 


'  .  _ _  .  *<2-34  ;  I 

-  -  -  - - —  ■■■' — ■ — - - - - - -  - : — ~-7r~T--~  -".rr,aUj; 


O  <" 

fi-  ' 


v-.«y  \  r-  • 

-  ■  V1 4A’ 

1  "  r&  ■  O  \  &  • 

o 

{-5V 


•  c  >  •'  4  ' -  ^  ^ 


LINE  »  0 

—  Write  tlb,45)  tpage 
WRITE  (LB, 44) 

24  WRITE  (LB  ,46)  RHO ,  AZ 1,  EL6V  *  RADAR  M  l,  - 

LINE  =  LINE  +  1 

25  CONTINUE  1  ,  _  !- 

IF(  JPR INT «LTa4  )  GO  TO  230 
IF(HIGH.GT.O. )  GO  TO  21 

NOR  BIT  »  NORBIT  -  1  "  ‘  Tw 

IF(  NORBIT  )  22*91*21  _  _  _ .  . 

”22  fNORBIT  »  NORBIT  +  2 

I F \  NORBITiEOvO  )  GO  TO  91  . . . 

READ  (LC)  PERIOD 
REWIND  LC 

SHIFT  =  SHIFT  ♦  PERIOD 
GO  TO  21 

48  FORMAT! 15, 4F1 5. 8, 15) 

' 43  FORMAT { 12F6.3 ) 

45  FORMAT!  72H 


FINA 


♦  *  40X*  6H  PAGE  *12) 

42  FORMAT ( I6,4F9*4) 

910  F0RMAT(13H0SUBSTAGE  N0.I2*17H  THRUST* INITIAL  *F10.2*14H  LBSjl 
1L  *F10. 2, 16H  LBS*,  FROM  T  «FT.3»11H  MIN  TO  T  -F7. 3  *4H  ('I N ) 

"911  FORMAT!  20H0  SPECIFIC  IMPULSE  *F7.4,14H  M2N.  NOZZLE  *F6.3,12H  SO* FT 
G  ~I2» 15H*  THRUST  ELEV  »F6,2,17H  DEG.*  AZ.CQRR.  «F7.2*6H  CLOCK) 

912  F0RMATM8H0  COEF.  OF  DRAG  *  F5i2*24H*  AREA  X)F  CRQSSECTIOH  «  FtU2*2 
+6H  SQ.TN.y  INITIAL  WEIGHT  «  F9iZ,21H  LBS*  EXTRA  SUBSTAGES  12) 

913  FORMAT ( 8H0MACH<  NO, 11F10.1/10H  CO  *li£l0c3)  „ _ 

601  FORMAT { 22H0INITI AL  POSITION  X  -  F9.3,  8H  *  V  *  ~F9.~3,’  8H  *'  Z  *  F9 

1.3*20H  NM,  INITIAL  TIME  «F8.3*4H  FUN) 

"602  FORMAT ( 24H0 INITIAL  VELOCITY,  VX  »  FL4*7*9H  *'  '  VY  *  F14.7*9H  ,  VZ 
1=  F14.7,7H  NM/MIN) 

961  FORMAT! 30H0INITIAL  GEODETIC  COORDINATES  FTV3V8H  NORTH*  F8.3*7H  EAS 
+T*  F7.3,27H;NM  HIGH,  INITIAL  TIME  *F8.3il6H  MIN#;  AXE_S  FIXE0I2)_ 
902  F0RKAT(29H0INITIAL  AIRSPEED,  DIR. ANGLE  F7<,2*i3H* INCLINATION  F6.2* 

+  11H  MAGNITUDE  F9.4,7H  NM/MIN  ,26X10HAXES  FIXED12) 

'  44  -FORMAT  <*1-20H  TIME  ‘  GLAT  ‘  &LdN*~““  HEIGHT*  .HEADING 

+  INCLINE  SPEED  RANGE  AZIMUTH  .  ELEY  SENSOR//) 

46  FORMAT!  75X#;3Fll.  3,  4XF6.G1 

47  FORMAT!  F9, 3, 9Fli,3*4XF6*0)- 

'  49  FORMAT ( /  9H  MISSILE  ,I4*11H  LANDED  AfYFY.3,T3H  DEliT^DRTH,  F87J* 

+  17H  DEG.  EAST  AFTER  F7.3,16H  MSN.  OF  FLIGHT.  I 
1015  FORMAT (23H00ID  NOT  COMPLETE  ORBIT)  “  ” 

4001  FORMAT!  1H  ,9HCR0SSING  I5»6X*5HTIME*E12.6*6X*3HL0HGfJ:_8_._3_)_ 

2016  FORMAT!  15H0N0DAL  PERIOD-  E12.6J  — . 

606  FORMAT ( 27H0INTEGRATI0N  STOPPED  AT  T  »F?.3*4H  MIN) 

•END'  '  "  ’  . 


SIBFTC  CMHlpy  DECK 

subroutine  cohand 

,  CQMMQN/PAft'/L  A  <Lfl  «PAT (4) iCQNV«SC&t AT iSCALONiPOt.  AT«PfllON  .. _  ' 

COMMON /FA C/NF.iNaMTYP  (30)  *MF£O0)  »DATT(60)  *¥11000)  «X  11060} 

_ _ COMMON /WHP/NW  «MWZ  ( 1  ?)  «  VW ( PQOO  )  f XW  (POOP)  ■  ,  ,  . . 

“  COMMON /MA P /NAME « CHAR i YN«YS«YSC 4 XV«XE»XSC*AMAP( 50*120 V  “  '  r 

DIMENSION  MWTYP(12J  .  ...  ,,  .  .. 

DATA  MWTYP  (1)  /72H  I  2  3  4  “g  <T~  T~  8  “ 

1  9  10  11  12/i  BLNK/6H  /♦ .  KW/iOOO/aKF/lO 


200 MY /SO  /  4MX/120/*  W/l.HW/  »F/1HF  /  ; 

READ  (LA«40)  YLAT<Xt  ON*nLAT«fltdN*Sl AT  *Si.flN 


-  -  -  •  .  *  u  ■  1  ”  ■  ■  TUI  vn  7  nwri  ■  T  ■  y  vn  T  ^  i  T 

'WRI?TE<L.?T.40)  YLAt»XLON»OUT4ttLON»StAt»SlON 
YN  s  srAlAT<>-(;l.»e0NV)OAtMTrrLATi  *  POUT 


*V  =  SCALON»(l«-C0NV)«AiNf  (XI.CiN)  «”A0t5N" 
YD  a  SCALATo(l.«CONV)«Ai^T(DLAfi 
XD  =  SC  At  ON# ( 1*-C0NV) #AlNT (dLON)  *  “ 

YSC  s  .SCALAT,#<SLAT-CONV#A!NT(SUT)l..  .  . 


XSC  =  SCAL0N#<SL0N-C0NV«A!NT,(SL0NU 
YS  s  YM-YD  . 1  .  ‘ 


XE  s  XV-XD 
NY  =  (Yn/Yscv+i; 


nx  =  (xn/xso  *1  i "  " 

I F ( N Y . 0 T M Y ; OR i N X i 6,T i MX)  RETURN 
DO  ’ 1  !Y=liNY 
DO  11  TX  =  iiMX 


TTTmAP  ( IY,!*)  =  8LNK 
12  READ  (LA* 4 1 )  NAMEVTY.CHAR 


WRITE! L'O  *A  1  >  N AM E  «  TY  ♦  CH  AR 
IFUY.NF.WL  GO  TO  13 


Call  compUt(nw,mwz «nwtyp»kw*yw»xw)v 

GO. TO  IP 


TXTF  (TY.NPiF)  GO  TO  14 

CALL  C0MPUT(NF,MFZ»NAMTVP.KF*Y*X> 


14 


GO  TO  12 

YMN  a  ATNT4YLAT) 

T  Y  M5~=  YMN"AINT (dlayT 
XMW  ,=  ATNT(XLON) 


Thz  =  XMW^AINT (nLON) ”  .  ~  : 

WRITFflP.,4?)  YMS»YMN,XMf»XMW  . . 

WRITE  (LP*43)  I  <  AMAP  UY*  IX)  *  tXPl'MX)  ilY»i«Nvr 

RETURN.  .  ....  .  ’  .  ’ 


xv 


40  FORMAT (PH  CARDLBL  6F8, 

41  F0RMAT(1XA6»1X2AU  .  -  v 

TTVTOP*rniniOT3i^ 

.1$  LATITUDE,  AND  F4»0t9H  AMD  F^OtlON  DEGREES  lONstTUftPrl 


FUNCTIONS 
«  ■  ■  . . . 

1.  RANDOM 

■  4  > 

2.  PKAY 

3.  DIST 

'  c 

4*  LOCATE 
5.  VIND/DPNV 

SUBROUTINES 

1.  RANORM 

2.  COORD 

3.  COMPUT 

4.  DUMMY 


'  "  ■*'  -  '  , 

,  '  6  *  4 

'  _  ...  , 


KraoreHSijf 


SIBFTC  RNn.Vn*  deck"  - - - - 

FUNCTION  RANDOM  (X) 

Common  /par/  datTaTTra 

no  1 9  J  a  1 .  p 


.  BANnp  =  RANDB  ♦  RANHA 
RANOA  a  RANT)B  -  RANDA 
T;F < R A N n 0 « G E *  1 . 0 )  R A N DB  a  RANDB  -  liO 
12  CONTINUE  ,  , 

Random's  randa  r  1  1 — - - — 

RETURN1 


:,;#V '  - 


The  PurP°se  of  the  function  RANDOM(X)  is  to  generate  a  I  "  :  -  *  • 

5  •  •'  ».  v  ■"]$ 

new  random  number  from  a  uniform  distribution  on  the  interval  0  to  1,.  V  ‘ 

«*  -  !  **  ’*  '•  •  ,*  <«•  A * 

The  procedure  used  is  a  Fabonachi  type  in  which  the  last  two  numbers,  'V  *\'  <' ^ 
generated,  RANDA  and  RANDB,  are  remembered  and. used  inform  ••  V  -  v  .1  *  1 

/  '  '  '  ,v  ^  \  .  : 

*  *  '  *  ?  o  ^  f  'T*> 

the  next  one.  Two  non-commensurated  numbers  are  supplied  frdim  ? 

/*  "  ^  ’  1.  '  '  4*'C 

the  compiled  parameter  list  to  start  the  process.  ~  v  .  ’  -  ^  tc 


BFTC  PK A?DK,  .  DECK _ _ 

^UNCTION  PK  AY (AT»AW*5Zi  '  .  * 

COMMON  /P Ap/  nAT(14) <CPKl«CPKgtCPK3  — — 

TTfMFNSlON  NC  (50)  _  .  . 

n A T A  ( Nr  t  J ) .  Js  1 « SO )  / 1 2*6 1  i S3 «»3 t37 « 


ARG'a  AW/ AT 

ANC  a  SZ/AT  ♦  CPK3 


ANC  =  AMAXl ( le  «ARG)  /ANC. 
PS  s  EXP(-ANC) 


100.«SGRT(ARG) 


NC  s  NC<K+ll  '  „  i  .  , 

(i^Pfi)»AHfNUl«eARO«CPta»ANC»»CPKg) 


PS*#M«/AN 
RETURN  . 


The  purpose  of  this  function  PKAY(AT,AW, SZ)  is  to  compute.  .  i  V\ 

’  *  “  '  .  '  o'/  *<?  '.''v 

-  , .  ^  ,,  •<  . 

an  estimate  of  the  average  fraction  of  remaining  value  destroyed  :by  ,  ; 

each  weapon  of  lethal  area  AW  and  variance  SZ  that. may  be  aimed  at  i\ ,  . 

«  ' , 

a  target  of  area  AT.  The  shape  factors,  CPK1,2,3  are  normally  get  „ 

- '  '  '  v 

so  that  for  targets  of  area  less  than  the  weapon  (aimpoints  all  at.  center  .  ^ 

'  ;  '  •  KM": 

of  target)  the  function  reduces  to  the  usual  single  shot  probability i :  /  Vi 


•  -  i  %  f* 

'  *  A- 

,  , ' 

.  V 

.  •?  v  ■  - 

Ll.  J.ii: 

'"■(f 

'  *  ~  VV‘ 

~r  •  u.~ 

.V  v 

:ert. 

& 

■A 

\  * 

.  '  Vi 

-  \ 

,.k  "  A-' 

<< 

;.-v  v.'  iy, 

*  , 

•'iv; 

other  shape  factors  i 


SI© FT C  DIS1DK  QCC<  ' 

FUNCTION  OTST(VA<XA«YB«XB) 

CoMmon  /par/  la.lb»rAd1us  • J "  "  •  -  ' ?  "  “ 

YAR  =  YA/aAOlUS  ' 

YSR  *  VB/RADJUS  "  “““  ‘  1  “ 

YR  =  (XA«Xft)/RAmUS  _ _ _ 

.COSANG  =  SIN  (YAR)  °SlN(YARl  ♦  eOS(VAR)«COS<Y®^)«COSlXRi 

nisr  a  ?, rraqiusoatan i sort ( ( 1  ♦■cosang) / Mc+cosangmi 

"  ft E til » N  .  '  '  _ 

FND 


The  function  DIST (YA ,  XA ,  Y B ,  XB)  computes  the  great  circle 
distance  in  nautical  miles  between  any  two  points  A  and  B  on  the  surface, 
of  the  globe.  The  latitude  and  longitude  of  each  of  the  points,  expressed 
also  in  nautical  .miles  (normal' machine)  form,  are  the  arguments. 
RADIUS  refers  to  the  radius  of  the  earth. 


im flh 


man 


S' 


-f 

■*;■•■ 


■;M 


;■  W.i 


SIBFTC  LOCm  !)FCK 

FUNCTION  LOCATE (V«X«K) 


y  *  ,  ■ 

common  /PAR/  LA,LB«RAD!US 

common  /FAC/  NF,NAMTYp<3Oi«MFZf3b)*VALT(3O!«&rtDTf30S* 

:  ■:  . 

*  uaooo)  ♦vuo'boi  «wuooo)  »R(ioooi«s(i6oo)*.Tfio6o) 

KR  s  1 

.1  )  KB  S  MF7CK-1)  ♦ 
MF7(K) 


=  y/ rad jus 
'  s'  x/radjOs 


=  COS (A j 
3  C0S(B)#THE0 


SHED  =  STN1B)  <*TH£0 

THED.  s  STM  (A)  -*  .  .  . 


DPL  =  CHED<»U(KE)  ♦  SHED«V(KEr ♦  THEDOWIKEJ 
DO  10  J  s  KR«KE 


TF {  DPL.GE.Oi)  GO  TO  lfl  , 


DPR.LT.O.)  GO  TO  10 
, DPM  =  CMEn»R(J)  „<>  SHeD*SCJ>  ♦  THED»TCJ>  • 
C  TO  ?3 


ao  locate 


ETIIRN 
END 


KB  +  1 


5  .  -  -  i 

r  4  • 

■  *■’  |*0  O'  0 

i !  •-  "i 


The  functionXQCATE(Y,X,k)  determines  in  which  section  of  *Va  ^ 

,  *  *  4  :  o'  ", 

a  mesh  a  particular  point  lies.  The  arguments  are  the  latitude:  and  ♦  : 
longitude  of  the  point  and  the  facility  type  number  K  of  the  desired  ^  4  -f 

meshV(usually  preselected  by  the  calling- program).  The  procedure  tests  .. 

•  V  '  *  f 

-  v  ✓  .  t 

all  sectors  in  order  from  one  to  last  until  a  positive  association  is 
determined'.  If  no  association  is  found,  the  point  is  automatically 
classified  as  last  one.  A  point  is  in  a  sector  when  it  lies  to  the  right  < 


V  if  •  •  *, 

'{!  .  V 

•  1  o 

?  *•"'  ft  5 

■f  , 


ueicnumeu.  xi  iiu  «sww«ou  is  luuiiu,  uic  pouic  is  auwmducauy  .-\  •>  '■j.fri* 

•  _  -  H. 

i  .  .  .  ’  -v 

classified  as  last  one.  A  point  is  in  a  sector  when  it  lies  to  the  right  '  ; 

*  ‘  '■  "  i 

of  the  left  hand  heading  looking  outward,  on  or  to  the  left  of  the  right.hand;(  ;  ^ 

5  /  ‘  -  '*-*$''*  ,J  ,A  "  "V"* 

heading,  and  on  or  exterior  of  the  section  arc y.  Exterior  means  to  the  ,  :  ! 

'  ,,  T  r‘r  ‘•i** 

^  lr  ,  ^  -ir-  „  tWi..  •.J'x.’’  — ■»]  *i  SL 

left  as  the  sector  posts  are  traversed  ih  increasing  order »  3*.^ 


A.  f  s  '  *  o' 

vV  * 

s>  O  V' 


' '"W' :i'P' 

•  *’'*  ’  -  v  ^ 

J*-  •  A 

>  ...  ^  4  I  -  1 

V” 


*  .  -  ;T-:  i,  '<•  •  -  -  •  <  . 

-  . ;  -  . 4  '  ..v-'-'-'' 

«-  o-  .  ,  '■ ' > C  «  A  ^  i 

-  •' ;  >;  .  .  ,•  .n 

*  '*  .  V,  *  *<#.*/&*&*> 


FUNCTION  VINDIX ,XQ) 
DIMENS  i  onT(T61\C‘(41 
DOUBLE  PRECISION  OP 


DATA  B/lV9i333333333,J5t.l.5,2.291&666Gj&,-2.45033333«Vl* 54166 
73 737 r7T237~17T2 5 ^  .375 ^I7^r6 666661^-1  *333333358 41  666'565 V 
+)  -*.833333333, *416666667/ 

~  TFT  XD.NE.0.:i  'GO  TO  70  — r— -  •  •  -  -  - -  7  : 

H  fc  0  ..... 

H"b  0  ‘  ‘  T 

TT  *  0  - 

do  i  j-=riT50“  ~  ~  :  --  . 

1  XX  (J).  k  0 

GO  .TO  120  -  “  “  ~ 

70  IFt  XO.EQ.H  .OR.  MVGTJ4  1  GO  TO  I74t6l»68,67|.6§^6S,63«.62).»M 

H  -~XD  :  ;  "  r‘~r!  ; 

M  fc  6 

cm  =*  btit  “  r  '  •  ••  .  ..//• " 

C  (2).  =  0  •  .  ‘ 

- - CC3f^-% - - - - - ~ - - — - “T— 

C<4)  =.  0  •  '  .  •  ■ 

~SO  TO  73~  ~  :  '  :  :  ■  ,  -  '  -  ‘ 

6S  ft  h.  M-i  .  •  .  ..... 

C(l)  =“BTM=2D  — - - - - - - - - — 

C  ( 2 )  --B  ( M-*2 ) 

GO  TO  120  •"  -  —  ■  — 

6 1  Mv  =.  7  '  -■  : 

C(  1)  *  TO)  '  “  ~  ~  :  1  :  “  ■“ 

C (2)  =  -B(3J  . ■  - 

DO  3d*  1/50  ~  . -  ■  "J'’-'-.' 

3  XDllJl  A  XD3 IJ)  .... 


GO  TO  74 


*  ,'t 


7V  ~  'r  T-l 


>V  A'f;, 


DO  4  J 


B<j4i3); 


N  =  8 


0  T  O  .73 
2 


GO  TO  145 


C(U)  B< J*4 ) 


,  v» ,  •*  -  4 

.»  "•* '  v7»  *.  /a  -iv^  v>  .of.  s'#* ?j 

.  1  O.  'f‘  *  W  «V*  rr.  ».  if  * 

\  ^ 
V‘  J*  $?*•"■  ',  <*AV  V  ,V  ' 

' .  T  ^ .  ’  *»  ,  <■«  «  *,  v 

« *  « V*  ^  * v 

"  f  ~'  •*  **  *  'K 

'  ; 


‘  ..  *':7.‘  ' 


103  I  -  1  :  - - * 

...  '  GO  TO  54 

ENTRY  O'PNVt a>XOJ"  -  : 

_ IF.l  I.GT.50  )  STOP 

R  *»  C(  4)  *XD3(  1)  “  ' 

^ _ :IF_t  K.EQ.1  );  X031II  «*  XD2ti). 

R:  m  R  *  CC3J*XP2<  fT“ 

IRC  K.EQ.l  );  X02(il),a«  XOim 

R  R  *  C  (2  ).*XD1  ( 1 1 
■  -■  IF L  K.EQ.l  );  XOlCfl)  n  XD 

R  te  ft  J-  C(i).tXO  r-r-r— 

opRm  *  x 

_  0PR(2>).  «»  xxtu;  “  ‘T~“ 

_ ..DP  *  (R*H)  *  OP 

xxm,  »  OPRfc2);  1  “  ~ 

_ -  i  >  i  +:i .  ■ 

54  VINO  =T  DPRU)  '  ■  

-  „ RETURN 

~  END  :  '  ~ ^ - - — " 


'  o"*  '.  •' 

, '  .  ,1-  -1  .  <  • 
;  '  o  £  o .. 


The  function  VINDfX.XD)  with  alternate  entry  point  *  7 

*  *  *  • 
DPNV{X,  XL')  is  a  generat;purpose  routine  to  integrate  up  to  fifty 

simultaneous  first  order  differential  equations .  It  Is  used  by  SURVEY 


V  *  *  .  ?  K  1 

,.v  *  • 


t''  '»  K 

r  -  * i  *  • 


differential  equations  governing  a  missile  or  satellite  orbit.  While 


'C,  .  t,  v  , 

,  ‘*"*H’"**U 


,  it  achieves 


,  then  retrelving  next  call 


the  second  word  worth  of  digits  for  each  variable,  A  call  to  VIND  Is 

*  » „  » 

used  to  increase  the  independent  variable  where  X  is  the  current  time 

i  *  t  '  T  ~ 

and  XD  the  increment!  a  call  to  DPKV  Is  made  for  each  dependent 


-  4 


*-t*»  »V«  *  , 


■  >?  ;  ;  ’  *  ;  ‘.-T  .'r,T 

, .  ‘  s'  *,l  ■  •  /  ,  ,  >  r\  I 

.X...  ,,^r,  .  ^  ‘  rv  '  ,'t’,  ■*. 

f  *  ,*  ^  '  '  t.  ’  *•'**  V*  *  /  { 

S  «  i  ■  v  *  '  ?  <  .  ‘  » 

'  '  'ff  *  *  1  '  *  11  -  /*  *4*  f  *Vl'  } 


i»£‘  t  > 


•  .  *  ■  .  .  1  •  '  -  „  ,  —  - 

;*  ■-  vv-.,  ,  -l  *  ; 


, ,  *  «►  t  -,v- 


\  •  o! 


,  •-**  .  -  V  .1-7.-.  i‘v  ItW,-.  .  •  *  v/  •  tl 

i  ‘  **,V  -  ^  '/  .  i  ^  J  '  Af  V  Jt**  *  >  '  •>  »  *  ,  V  /■)  j  .it 

. -  '  >%'v '  f /  -V  v*'  '•  H  I| 


•  ,  v!,-w 

■*!/-'  i .  *  ,  iXjiL  4 


.  V>  ^  -  -V  *,<  v  n,  9r  v'i  V  VA  V'*  A* 


,  ,  -  fv>-  ~-  v  r,2  -,-<*»■  r  „ 

*  >  O 


x>  -  I  •'  ■  ..  -•'-.*  .'■_  r>--  ”  ''•'  -  \  ,^X-.'  '  «!'7\<  t„?6  3 


4.-0, 
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'o  <,  •ii,;A  ■";  ‘-<',  v*1 
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■  -  ,  ■•*  '-v^f  - >  •/.<>' '; 

-  "  „  ct  ;  '«a t'y'  '-.  vV 

.  «?  o  *  >  >  -  -  $  ..***  \<s** 

/•  <* 

variable  to  be  integrated  where  X,  XD  are  the  current  values  andiirst 


.  *?  0  *’  \  '*■'  *< 

o'-  *  [  \  >n  -  /  *  «r<< 

'  v  b 


derivative  respectively. 


,  .*  ,*vis 


*  V*- -  .*-■ 

,?  £  .-.'V'.’Vi 


This  program  was  originally  obtained  as  a  biliary  deck  in  70 4y-  ‘  v-'l 

*  •  v  -X  •  V.4.,  '  X  S^l 

,  ,  ..  f  »  ?X«S«'  %'  ,  lX  •, 

machine  language  from  Lincoln  Laboratories  for  use  with  Fortran  II  o  '  ■  0/  .  •  -  ;, 
programs.  It  was  rewritten  by  the  author  in  MAP  language  to  be  -  ,  ..  V  ^••vv;‘r^ 


compatible  With  the  Fortran  IV  computer.  It  is  this  form  or  its 


.  c  x  <>;<.  .  •  v‘*,b  *  ? 

C  >  r  <v  ’  .\j}  ,  .  -v 


'W'o-V 


-.-n-'.y.. 


'V.  V 


in  this  report,  however  MAP  programs  c. e  hot  easy  to  understand. 


Sin-e  the  MAP  version  has  been  rewritten  in  Fortran  IV  under  another  *  -  ^ 

contract  (to  make  it  useable  on  a  GDC  3600. machine)  this  more  readable 
version  is  used  here.  It  has  been  tested  on  the  7094  with. SURVEY. 
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